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1  Introduction 


A  companion  [1]  to  this  Report  discusses  motivation  for  the  analysis  of  near  net-shape  stress-coated  mem¬ 
branes,  and  presents  details  of  the  derivations  of  four  different  theories  of  the  mechanics  of  such  a  laminate 
using  the  method  of  asymptotic  expansions.  Our  goals  here  are  to  present  solutions  of  selected  boundary 
value  problems,  and  to  compare  the  results  with  those  obtained  by  the  finite  element  (FE)  method  where 
available.  Each  of  the  next  four  Sections  will  be  organized  in  a  similar  way,  i.e.,  the  governing  equations  for  a 
given  theory  will  first  be  restated  from  Reference  [1],  boundary  conditions  will  then  be  imposed,  and  details 
of  the  solution  procedure  for  each  such  set  of  boundary  conditions  will  be  presented.  Finally,  graphical 
comparisons  will  be  made  between  the  solutions  of  a  given  boundary  value  problem,  and  those  using  FE 
analysis. 

With  regard  to  nomenclature  and  notation,  all  problems  presented  here  involve  an  JV-layer  structure 
consisting  of  a  membrane  substrate  of  thickness  hx  —  hs  having  N  ~  1  coatings  with  thicknesses  hi,  i  = 
1, . . . ,  JV  —  1,  and  total  thickness  h  =  hs  +  hc,  where  hc  =  J2k=i  hk  is  the  total  coating  thickness.  This 
structure  has  a  circular  boundary  of  radius  a.  The  reference  configuration  S  of  the  coated  membrane  is 
defined  by  a  mapping  from  a  purely  geometrical  region  C,  which  we  refer  to  as  the  reference  placement.  This 
mapping  is  illustrated  for  a  single  coating  in  Figure  2  of  Reference  [1],  and  reproduced  below  as  Figure  1. 
As  shown  in  the  Figure,  the  coating  and  membrane  thicknesses  are  assumed  to  be  constant  along  any  line 
parallel  to  the  axis.  The  reference  placement  for  N  —  1  coatings  is  illustrated  in  Figure  5  of  Reference  [1], 
and  reproduced  below  as  Figure  2.  Note  that  the  coatings  are  “on  the  bottom”.  The  basis  of  this  rather 
counterintuitive  orientation  is  the  optical  convention  that  the  direction  of  propagation  of  light,  which  strikes 
the  coating  first,  coincides  with  the  positive  Z-direction  (upward  in  our  Figures).  The  coordinate  £  =  Z+h/2 
shown  in  Figure  2  is  convenient  for  labeling  the  top  and  bottom  surfaces  of  each  layer,  and  we  note  from 
Figure  2  the  important  relation  £*  =  Y?k= i  hk- 

The  following  constants  appear  in  the  through-the-thickness  integrals  of  the  constitutive  relations  defining 
the  stress  resultants  and  stress  couples,  treated  in  detail  in  Appendix  A  of  [1]: 


N  H  N- 1  N 

v  =  e hiSi ’  M  =  2  £  E  h -  Sk )«  (i-i) 

i~l  i—  1  k  =  i+ 1 

N  N 

A  =  Y,  hiQi,  Ay  =  E  (1-2) 

t=l  i= 1 

N- 1  N  .  N-l  N 

B  =  ^Ti  E  hihk(Qi-Qk),  £„  =  -  £  E  hihkiQit'i-Qki'k),  (1.3) 

z  —  1  k  =  i+1  '  *=1  k  —  i+1 

N  -  N-l  N 

A&  =  E  h>Gi’  5e  =  -^  E  hihk  ^Gi  ~  Gfc)’ 

t=l  i=l  k  =  i+1 

1  N  r  , 

D  =  12  E  &  hilh2i  +  3  -  /li)2  -  12^-1  ( h  -  6)J  ,  (1.5) 

i= 1 
1  N 

Dv  =  jzY'QiUt  hi  Ih*  +  3  (h-  hi)2  -  12&_i  (h  -  €*)  ,  (1-6) 

Z  i=l 
1  N 

=  Jz  E  Gi  hi  [h*  +  3  (b  -  -  12^-1  ( h  -  6)]  -  (1.7) 

i= 1 


i 


Surface  of  Mold 


Figure  1:  Definition  of  the  reference  configuration  S  (upper  part  of  Figure)  of  a  coated  membrane  shell  of 
revolution  as  a  mapping  from  the  reference  placement  C  (lower  part  of  Figure),  assuming  the  layer  thicknesses 
to  be  constant  along  any  line  parallel  to  the  axis. 


where,  in  (1.1),  S,  is  the  residual  stress  load  in  layer  i.  The  elastic  constants  of  the  material  in  each  layer 
appear  in  the  forms 


Qi  = 


Ei 


Gi 


Ej 

1  +  Vt’ 


(1.8) 


where  Ex  and  vt  are  Young’s  modulus  and  Poisson’s  ratio  for  material  i.  The  areal  density  of  the  coated 
membrane  (mass  per  unit  area  of  the  circular  disk  perpendicular  to  the  axis;  units  of  kg/m2 )  is  defined  and 
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Figure  2:  Reference  placement  C  of  an  iV-layer  stack. 


denoted  by 


N 

7o  =  ^2  hipoi, 
i= 1 


(1.9) 


where  poi  is  the  volume  mass  density  of  the  material  in  layer  i.  The  external  loads  are  the  pressure  difference 
p  =  p~  -  p+  between  the  two  faces  of  the  coated  membrane,  and  the  gravitational  body  force  70  g  (where  g 
is  the  gravitational  acceleration)  in  the  positive  axial  direction  (up)  of  our  two  Figures. 

We  realized  only  after  the  publication  of  Ref.  [1]  that  the  nine  A,  B ,  and  D  coefficients  were  not  inde¬ 
pendent.  We  in  fact  find  the  following  three  relations  between  the  multilayer  coefficients  (recall  that  each 
layer  is  assumed  to  be  uniform,  homogeneous,  and  isotropic  in  its  material  properties): 


A@  =  A  -  A„,  Be  =  B  -  Bu,  D©  =  D  -  Du. 


(1.10) 


Similar  to  Jones  [2,  p.  155],  we  refer  to  A  as  the  extensional  stiffness ,  B  as  the  coupling  stiffness ,  and  D 
as  the  bending  stiffness  of  the  multilayer  laminate.  It  is  convenient  to  also  introduce  extensional,  coupling, 
and  bending  Poisson’s  ratios  defined  by 


Av  Bv  ,  Dy 

va  =  -j,  vb  =  -g~,  and  vD  = 


(1.11) 
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respectively.  Note  that  for  a  single  layer  with  Poisson’s  ratio  u,  or  for  a  multilayer  in  which  all  layers  have 
the  same  Poisson’s  ratio  v  (an  unlikely  possibility),  the  three  Poisson’s  ratios  defined  in  (1.11)  are  all  equal 
t°  v  as  well.  Equations  (1.10)  can  be  written  in  terms  of  the  stiffnesses  and  Poisson’s  ratios  as 

Ae  =  A  (1  -  i/A) ,  Be  =  B  (1  -  vB) ,  De  =  D  (1  -  vD) .  (1.12) 

The  most  important  initial  shape  for  optical  purposes  is  a  paraboloid,  so  throughout  this  Report  we 
consider  a  membrane  cast  on  a  paraboloidal  mandrel  (also  referred  to  as  “the  mold”),  then  coated  on  the 
mandrel,  and  subsequently  released  from  the  mandrel  as  a  near  net-shape  coated  membrane  laminate.  We 
assume  that  the  coated  membrane  middle  surface  is  governed  by  the  equation 

m  =  jj  («>  -  *!) ,  (M3) 

where  /  is  the  focal  length  of  the  reference  paraboloid.  A  flat  middle  surface  is  modeled  by  a  paraboloid 
having  an  infinite  focal  length,  in  which  case  T(R)  =  0.  From  (1.13),  the  central  (or  vertex,  or  apex) 
displacement  of  the  paraboloid  is  given  by 


(1.14) 

(1.15) 

(1.16) 


r0  _  1 

a  8 F*  ‘  (U7) 

In  the  optics  literature  a  paraboloid  with  (for  example)  an  /-number  of  2  is  referred  to  as  an  //2  paraboloid, 
and  we  indicate  this  by  writing  F*  =  2.  The  ratio  (1.17)  for  an  // 2  paraboloid  is  1/16,  or  0.0625. 

It  is  often  more  convenient  to  write  (1.13)  in  terms  of  a  parameter  k  defined  by 


1  1 

“  2/  “  4 aF#' 


(1.18) 


in  which  case  we  have 

r(R)  =  f  (°2  -  r2)  •  (1.19) 

Note  that  the  principal  curvatures  m  and  /c2  of  the  paraboloid  are  given  by 


M  (  1  -  K2R2  )3/2  ’  aild  K2  ~  (1  _  «2iJ2)l/2>  C1-20) 

which  are  equal  only  at  the  vertex  R  =  0  of  the  paraboloid,  where  m  =  k2  =  -k.  Thus,  «  is  the  (negative) 
vertex  curvature  of  the  paraboloid,  and  k  — ►  0  as  /  — t  oo,  i.e.,  as  the  paraboloid  flattens  into  a  plane. 

Unless  otherwise  stated,  all  variables  appearing  in  this  work  are  leading  order  functions  of  the  cylindrical 
coordinates  (R,  0,  Z)  on  the  reference  placement  C.  The  subscript  (0)  used  to  denote  leading  order  variables 
in  Reference  [1]  will  be  omitted  in  all  that  follows.  Reference  [1]  should  be  consulted  for  further  details  on 
the  definitions  of  the  functions  appearing  in  the  equations  of  this  Report. 
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2  Geometrically  Linear  Membrane  Laminate  Problems 

The  simplest  of  the  theories  presented  in  Reference  [1]  is  that  of  a  geometrically  linear  coated  membrane 
laminate.  It  is  the  only  theory  for  which  we  consider  solutions  that  are  non-axisymmetric.  The  equations 
governing  the  mechanics  of  such  a  laminate  are  given  in  §11  of  that  Report,  which  are  repeated  here,  beginning 
with  the  strain-displacement  relations,  then  the  integrated  through-the-thickness  forms  of  the  constitutive 
relations,  and  finally  the  equilibrium  equations  stated  in  terms  of  stress  resultants  (stress  couples  do  not 
appear  in  a  membrane  theory).  Throughout  this  Report,  a  reference  to  an  equation  in  [1]  will  be  distinguished 
by  appending  a  “.I”  extension  to  the  referenced  equation  number. 

2.1  Pressurized  Membrane  Laminate  with  ©-Dependent  Boundary 

For  membrane  laminates,  the  leading  order  displacement  components  are  given  by 

UR  =  u,  Ue=  v ,  Uz  =  w,  (2-1) 

where  u,  v,  and  w  are  functions  of  R  and 
in  equations  (11.20.1)— (1 1.22.1) ,  viz., 

4© 


The  stress  resultants  (11.29.1)— (11.31.1) 
constitutive  relations  (11. 17. 1) -(11. 19. 1),  yielding 

Nr  —  Af  4-  Ac*rr  4-  Av  e©©,  (2.5) 

Nrq  =  A@c°r@ ,  (2.6) 

Ne  =  Af  4-  Aue°RR  4-  A £00.  (2.7) 

Finally,  recalling  the  definition  t o{R)  =  T(R)  +  w(i?,0),  the  leading  order  equilibrium  equations  are  given 
in  terms  of  the  stress  resultants  by  equations  (11. 26. 1) -(11. 28. 1): 

(RNr)  r  -  iV©  +  NRe,e  =  0,  (2.8) 

{R2Nr@)  r  +  RNe,e  =  0,  (2.9) 

j^R  (f,r  Nr  +  w,r  Nr  4-  R  (f,rNrq  4-  w^rNrq  4-  --jfNe)  ^  +  {p  +  ^g)R  =  0.(2.10) 

We  seek  solutions  of  (2.2)-(2.10)  for  an  initally  flat  coated  membrane,  so  that  T^r  =  0,  and  with  u  =  v  =  0 
for  all  R.  When  these  hold,  we  have  e^e  =  e°RR  =  £00  =  0,  hence  Nr  =  iV©  =  Af  and  NR@  =  0.  Equations 
(2.8)  and  (2.9)  are  then  identically  satisfied,  and  (2.10)  reduces  to 

(RwtRAT)R  +  (^V)  q  +  (p  +  Joff)R  =  0.  (2.11) 


0  only.  The  leading  order  strain-displacement  relations  are  found 


If  u,e~v\ 

=  2  \v,fl  R  /  ’ 


(2.2) 


eRR  ~  U,R> 
v©  +  u 


,0  ... 

ee©  - 


R 


(2.3) 


(2.4) 


result  from  through-the-thickness  integrals  of  the  leading  order 


5 


Carrying  out  the  differentiations  in  (2.11),  we  obtain 


w,RR  + 


R 


R? 


(2.12) 


noting  that  Af  is  constant.  The  left-hand  side  of  this  equation  is  recognized  as  the  Laplacian  of  w  in  plane 
polar  coordinates,  hence  the  equation  can  be  written  as 


V2ie  =  - 


/o’ 


P  +  Jo  9 


where  /0  is  a  constant  with  units  of  length  (m). 

We  seek  solutions  of  (2.13)  in  the  form  of  a  Fourier  series: 


(2.13) 


OO 

w(R,Q)  =  [An(R)  cos  (nQ)  +  Bn(R)  sin  (nQ)],  (2.14) 

n=0 

where  the  functions  An(R)  and  Bn(R)  are  to  be  determined.  Taking  the  required  partial  derivatives  of  (2.14) 
and  substituting  the  results  in  the  left-hand  side  of  (2.13),  we  obtain 


{K(«)cos(n0)  +  B"(fi)sin(n0)]  +  i  [A'n{R)  cos  (nQ)  +  B'n{R)  sin  (n®)} 

n= 0  ^ 

71^  ^  ] 

“  -jp  [^n(-R)  cos  (nQ)  +  Bn(R)  sin  (n0)j  j  =  ,  (2.15) 

where  a  “prime”  denotes  a  derivative  with  respect  to  the  argument  of  a  function.  Since  the  trigonometric 
functions  sin  (n0)  and  cos  (nQ)  are  linearly  independent  for  each  value  of  n,  the  last  equation  yields  the 


following  ordinary  differential  equations  for  An(R)  and  Bn(R ): 

A'c!(R)  +  ±A'0(R)  =  -1  n  =  0,  (2.16) 

K(R)  +  |  K(R)  -  §  MR)  =  o,  ii  /  o,  (2.17) 

K(R)  +  |  K(R)  -  ^  Bn(R)  =  0,  n  /  0.  (2.18) 

Equation  (2.16)  is  easily  solved  to  obtain 

A0(R)  =  a0  -  ~R2,  (2.19) 


where  ao  is  an  arbitrary  integration  constant,  and  the  other  integration  constant  was  set  equal  to  zero  to 
insure  regularity  of  A0(R)  at  R  =  0.  Equations  (2.17)  and  (2.18)  are  “equidimensional”  [3,  §1.6]  linear 
differential  equations.  They  can  be  reduced  to  linear  differential  equations  with  constant  coefficients  by 
introducing  a  new  independent  variable  £  defined  by 

£  =  In  R  =>  R  =  e*.  (2.20) 

Setting  An(R)  =  j4n(e^)  =  ^4n(0  and  Bn(R)  =  Bn(e^)  =  Bn(£),  and  using  the  chain  rule  to  compute  the 
derivatives,  equation  (2.17)  reduces  to 

Kd)  -  n2I„(0  =  0,  n  /  0,  .  (2.21) 
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and  similarly  for  (2.18).  The  two  linearly  independent  solutions  of  (2.21)  are  An(0  =  en'  and  .?„(£)  =  e 
hence  An{R)  =  Rn  and  An(R)  =  R~n,  so  the  general  solutions  of  (2.17)  and  (2.18)  are  linear  combinations 
of  these,  viz., 

An(R)  =  anRn  +  cn  R~n,  and  Bn(R)  =  bnRn  +  dnR~n,  (2.22) 

where  on,  bn,  c„,  and  dn  are  arbitrary  constant  coefficients  depending  on  n.  However,  since  the  solution 
must  be  regular  at  R  =  0  for  any  n  /  0,  we  must  set  cn  =  dn  -  0  for  each  n.  Substituting  the  results  into 
(2.14)  yields  the  general  solution: 


1 

w(R,  ©)  =  a0  -  —  R2  +  [ anRn  cos  (n0)  +  bnRn  sin  (n©)] ,  (2.23) 

4fo 

where  the  coefficients  a0,  an  and  bn  must  be  determined  from  the  boundary  conditions.  Typically,  one  is 
given  the  form  of  the  displacement  field  on  the  boundary  edge  R  =  a,  i.e.,  w(a,Q)  =  W(@).  The  Fourier 
series  coefficients  C0=a0-  n2/(4/0),  C„  =  anan,  and  Dn  =  bnan  are  determined  in  the  usual  way  [3,  §5.11] 
from 


(2.24) 


w(a,  0)  =  W(0)  =  C0  +  ^2  [Cn  cos  (n©)  +  Dn  sin  (n0)] , 

n~  1 

Co  =  hCw{Q)de’ 

1  /*27r 

Cn  -  -  W(0)  cos  (n0)d0, 

ft  Jo 

1  f27r 

Dn  =  -  /  VF(0)  sin  (n0)d0. 

^  Jo 

Having  determined  Cn  and  Dn,  we  can  replace  an  =  Cn/an  and  bn  =  Dn/an  in  (2.23),  and  the  solution 
takes  the  final  form 

i  00 

w(R,Q)  =  Co  +  ~r  (a2  -  R 2)  +  Y  [Cn(R/a)n  cos  (n0)  +  Dn(R/a)n  sin  (n©)] .  (2.25) 

4/o 


n=l 


2.2  Vibrations  of  a  Membrane  Laminate 

We  conclude  this  brief  Section  with  a  digression  on  the  subject  of  membrane  dynamics.  We  begin  by 
asserting,  without  proof,  that  the  dynamical  leading  order  momentum  balance  equations  are  the  same  as 
the  equilibrium  equations  (2.8)-(2.10),  with  the  exception  that  the  right-hand  sides  are  replaced  by  inertia 
terms  comprised  of  products  of  the  areal  mass  density  70  and  an  acceleration  component  (each  equation  has 
also  been  premultiplied  by  R ),  i.e., 


(RNr)r  —  Nq  +  IVfle.e  =  70  Ru,tt, 


(2.26) 


(J?2lVft©)fi  +  RNe,@  =  loRv,tt, 


(2.27) 


[f?  (t,rNr  +  W'R  Nr  +  R  +  (r ,rNrq  +wtRNR@  +  “j^-^©)  0 

+  (p  +  'Tog)  R  =  70  Rw,tt-  (2.28) 
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We  consider  free  vibration  solutions  of  (2.26)-(2.28)  for  which  u  =  v  =  0,  assuming  there  is  no  pressure 
load,  and  ignoring  the  effects  of  gravity.  We  further  assume  that  we  have  an  initally  flat  reference  configu¬ 
ration,  so  that  r, .R  =  0.  The  conditions  u  =  v  =  0  again  imply  e°Re  =  e°RR  =  e^e  =  0,  hence  NR  =  Ne  -  Kr 
and  NRq  =  0.  Equations  (2.26)  and  (2.27)  are  then  identically  satisfied,  and  (2.28)  reduces  to 

(R wtRAf)R  +  JV-)  =  7o  Rw,tt-  (2.29) 

Assuming  Af  to  be  constant,  we  carry  out  the  differentiations  in  (2.29)  to  obtain 


The  left-hand  side  of  (2.30)  is  the  Laplacian  of  w,  hence  we  can  write  the  equation  as 


(2.30) 


V2w 


which  is  the  homogeneous  wave  equation  for  an  axial  displacement  having  propagation  velocity 


c  = 


(2.31) 


(2.32) 


The  solutions  of  (2.31)  are  well-known  (see,  for  example,  [4,  p.  635]),  and  will  not  be  repeated  here.  However, 
the  angular  frequencies  of  vibration  are  given  by 


_  c 

Umn  —  &mn  ~ 

a 


ftrnn 


(2.33) 


where  a  is  the  membrane  radius,  and  amn  is  the  mth  positive  zero  of  the  ordinary  Bessel  function  Jn  of 
order  n  (n  =  0, 1, 2, . . .),  i.e.,  amn  is  a  solution  of  the  equation  J„(amn)  =  0.  Note  that  the  tension  T  which 
appears  in  the  small  oscillation  theory  of  a  single  material  “drumhead”  is  replaced  in  this  theory  of  a  coated 
membrane  laminate  by  the  net  residual  stress  resultant  Af  defined  in  equation  (1.1). 


3  Geometrically  Nonlinear  Membrane  Laminate  Problems 


The  leading  order  strain-displacement  relations  are  given  by  equations  (9.27.I)-(9.29.I),  which  contain  non¬ 
linear  terms  in  the  axial  displacement  derivatives: 


n, e  -  v 
R 


+ 


T.rw  e 
R 


R  J’ 


(3.1) 


eRR  -  U,R  +  LftlUft  +  -  (tl), ft)2 


(3.2) 


The  stress  resultants  and  equilibrium 
as  those  of  the  previous  §2: 


-  M +  u  ,  Ke)2 

ee  R  +  2R2  ' 

equations  are  given  by  equations  (9.33.1)— (9.38.1) , 


(3.3) 


which  are  the  same 


Aft  —  AT  +  AeRR  + 


(3-4) 


8 


Nrq  =  A©  e°R&, 


(3.5) 


Nq  =  Af  +  AvCftj j  +  A  egg.  (3-6) 

and 

(RNr) >r  -  Ne  +  NR@,e  =  0,  (3.7) 


[R2Nrq)r  +  RNe,e  =  0, 


(3-8) 


[i?  (t,rNr  +  w^Nr  +  +  w,rNrq  +  ~^~N®)  Q  +  (P  +  7off)^  -  0.  (3.9) 


3.1  Reduction  to  an  Axisymmetric  System 

We  specialize  immediately  to  the  case  of  an  axisymmetric  system  by  assuming  that  none  of  the  variables 
depend  on  the  angular  coordinate  0.  Thus,  all  terms  involving  partial  derivatives  with  respect  to  0  vanish, 
leaving  the  following  simplified  set  of  strain-displacement  relations  and  equilibrium  equations: 


4e  -  \  E)  » 


eRR  =  U,R  +  F,RW,R  +  -  (w,R)  , 


(3.10) 

(3.11) 


(3.12) 


(RNr),r  -Ne  =  0, 


(3.13) 


(R2NrQ)'R  =  0,  (3.14) 

[i?  (T'RNr  +  wtR Nr )]fi  +  (p  +  70 g) R  =  0.  (3.15) 

Equations  (3.14)  and  (3.15)  can  be  integrated  once  to  obtain 

R2Nrq  =  Crq,  (3.16) 


and 

B? 

R(TtR  Nr  +  w,r  Nr)  4-  (p  +  70  9)  (3.17) 

where  Crq  and  Cr  are  arbitrary  integration  constants.  We  must  set  both  of  these  constants  equal  to  zero 
to  insure  that  both  Nrq  and  Nr  are  regular  at  the  origin  R  =  0.  Thus,  (3.16)  and  (3.17)  reduce  to 

Nrq  =  0,  for  all  R.  (3.18) 

R 

^,rNr  +  W'RNr  +  (p  +  jog)  —  =  0.  (3.19) 
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Prom  (3.18),  (3.5)  and  (3.10)  it  then  follows  that 


.0 

€h© 


(3.20) 


since  Aq  is  in  general  nonzero.  Equation  (3.20)  is  easily  solved  to  obtain  i?(jR)  =  vqR ,  where  vq  is  an  arbitrary 
constant.  We  observe  that  if  v  vanishes  for  any  non-zero  value  of  R,  say  Ro  ^  0,  then  v(Ro)  =  v0Rq  =  0, 
hence  Vo  =  0,  and  so  v  vanishes  for  all  R. 

The  system  of  equations  that  remains  to  be  solved  is  then 


eRR  -  u,R  +  r,RwtR  +  ^(w.r)2,  €ee  = 

(3.21) 

Nr  =  N  +  Ac°rr  -f  Au€%e, 

(3.22) 

Ne  =  N  +  Au^rR  +  Ae@0, 

(3.23) 

C RNr \r  -  Ne  =  0, 

(3.24) 

R 

r  ,rNr  +  w,rNr  +  (p  +  jog)  —  =  0. 

(3.25) 

3.2  Generalization  of  Hencky-Campbell  Theory  to  a  Membrane  Laminate 

We  continue  the  reduction  by  specializing  now  to  an  initially  flat  coated  membrane,  so  that  T  =  0,  reducing 
the  system  of  equations  to 


4a  =  «,a  +  \  i™,R? ,  4e  =  ^  (3.26) 

Nr  =  N  +  Ae°RR  +  Au€%q,  (3.27) 

Ne  =  AT  +  A„  e°RR  +  ,44©,  (3.28) 

(RNr),r  -  Ne  =  0  =>  Ne  =  Nr  +  RNr<r,  (3.29) 

R 

W,RNR  +  (p  +  Jog)  =0.  (3.30) 


This  system  is  augmented  by  the  usual  compatibility  condition  obtained  by  eliminating  u  between  the  two 
strain-displacement  relations  (3.26),  yielding  (see,  for  example,  [5]): 


Rcee,R  -  c°rr  ~  4e  -  2  • 


(3.31) 


We  proceed  by  introducing  a  dimensionless  coordinate  p  and  dimensionless  displacement  components  u 
and  w  defined  by 

R 

(3.32) 


P  = 

a 


^  u 

u  = 

a 


=  T. 
a 
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(3.33) 


together  with  the  following  dimensionless  constants: 


v A  = 


A  ’ 


E  =  f(  1-vl), 


.  Af 

T  ~  Eh’ 


Q  = 


( P  +  _7off)Q 

Eh 


and  two  new  dimensionless  dependent  variables  xr  and  xe  defined  by 


_  Nr-  Af 

Xr  —  —  j 

Eh 


_  Ne-AT 
Xe  ~  Eh 


(3.34) 


Note  that  va  is  the  extensional  Poisson’s  ratio  defined  earlier  in  (1.11),  while  E  acts  as  the  “effective” 
Young’s  modulus  of  the  coated  membrane  laminate.  Equations  (3.26)-(3.31)  can  be  rewritten  in  terms  of 
these  dimensionless  quantities  as 


1  u 

£rr  =  2,p  +  2  ™?p’  4e  =  -» 

(3.35) 

xr  =  _  v2  j  ieRR  +  VA  eee) ) 

(3.36) 

xe  =  Q_j/2)  (e©e  +  va£rr), 

(3.37) 

Xe  =  Xr  +  pXTtP  =  (pXr)iP, 

(3.38) 

W,P  (t  +  Xr)  +  q  |  =  o. 

(3.39) 

Peee,p  =  eRR  ~  eee  —  2^’ 

(3.40) 

These  equations  have  precisely  the  same  form  as  Campbell’s  [6]  equations,  which  are  themselves  modifications 
of  Hencky’s  [7]  membrane  equations  to  allow  for  an  initial  tension  in  the  pressurized  membrane.  Equations 
(3.36)  and  (3.37)  are  easily  inverted  to  obtain 


e°RR  =  Xr  -  VAXff,  (3.41) 

e|e  =  xe  -  VAXr,  (3.42) 


which,  together  with  (3.38),  can  be  substituted  in  (3.40)  to  get  the  compatibility  condition  in  terms  of  xr 
and  xe : 

p(xr  +  xe)>p=  -\w%.  (3.43) 


We  now  use  Equation  (3.39)  to  eliminate  w>p  in  equation  (3.43),  yielding  the  following  equation  involving 
only  the  x’s: 


<r + + + ?;  =  0. 

L  p  J  8 

(3.44) 

It  is  convenient  here  to  rescale  r,  xr.  x$,  and  w,  i.e.,  we  set 

T  =  ig 2/3  r,  xr  =  ig2/3xr,  xe  =  jq2/3xe, 

4  4  4 

w  =  q1^  w 

(3.45) 

to  write  equations  (3.39)  and  (3.44)  as 

wtP  (r  +  xr)  +  2p  =  0, 

(3.46) 
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and 


(r  +  Xr)2 

respectively. 

We  assume  a  power  series  solution  for  £r(p),  and  find  after  some  numerical  experimentation  that  only 
even  powers  of  p  will  contribute.  It  is  convenient  to  write  the  series  in  the  form 

1  +  £  ftnp2n)  •  (3-48) 

n=l  / 

From  this  and  the  second  equality  of  (3.38)  it  follows  that 


(xr  +  xe) 


+  8  =  0, 


(3.47) 


Xe  =  (px  T)p  =  -T  +  b0 

The  sum  of  the  last,  two  equations  yields 

xr  +  xe  =  -2?  +  2  b0 

hence,  noting  that  f  is  a  constant, 

(xr  +  xe) 


1  +  £  (2n  +  l )b2nP2n 


n=l 


1  +  ^  (n  +  1)  &2n  P2n 


n=  1 


=  4  b0  ^2  n  (n  +  1)  hn  P 


,2n-2 


n=l 


(3.49) 


(3.50) 


(3.51) 


This  expression,  together  with  (3.48),  is  now  substituted  in  (3.47),  and  the  coefficients  b2n ,  n  >  1,  are  then 
determined  by  equating  to  zero  the  coefficients  of  like  powers  of  p.  We  find,  using  the  Mathematica  computer 
algebra  system,  that  b2n  is  given  in  terms  of  b0  by 


i  _  ^2n  .  , 

■ 2n  ~  uZn  >  n  — 

°0 

where  the  /32n  are  purely  numerical  coefficients  given,  for  1  <  n  <  9,  by 

37 
27’ 

51523763 
5143824  ‘ 


&  1,  ft  -  3,  ft  -  ft  -  ft0 

_  1205  Q  _  219241  0  6634069 

567  ’  Pl 4  63504  ’  Pl6  ~  1143072’  Pl8 


(3.52) 


(3.53) 


In  order  to  determine  the  coefficient  bo  (in  terms  of  which  all  the  other  coefficients  i>2n  can  be  calculated), 
we  must  impose  boundary  conditions. 

We  consider  here  a  clamped  boundary,  requiring  u(a )  =  0  and  w(a)  =  0,  or  equivalently,  2(1)  =  0  and 
trj(l)  0.  From  the  second  equation  of  (3.35),  together  with  equation  (3.42),  we  have  the  following  power 
series  representation  for  2: 


u(p)  =  p4>e  =  p(xe  -  VAXT)  =  -q2/3p 


-t(1  -  vA)  +  ft  (1  ~  vA  )  +  b0  ^  (2  n  +  1  -  i/A  )  b2np2n 


n=l 


(3.54) 

Applying  to  this  expression  the  boundary  condition  on  2  at  p  ~  1  yields  an  equation  which  must  be  solved 
numerically  for  60‘ 


bo 


1  —  VA  +  ^  (2n  +  1  —  Va  )  &2n 


n~  1 


-  (1  -  1/A  )t  =  0. 


(3.55) 
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After  determining  60>  which  in  general  can  be  seen  from  (3.55)  to  depend  on  both  va  and  r,  the  series 
solution  for  the  radial  displacement  u(p)  is  given  by  (3.54). 

The  dimensionless  axial  component  of  displacement  w(p)  can  be  obtained  by  assuming  an  even  power 
series  of  the  form 

oo 

w(p)  =  E  C2^2”’  (3-56) 

n= 0 

the  derivative  of  which  is 

oo 

«>,p  =  2'£nc2np2n~1.  (3.57) 

n=0 

This  is  substituted,  together  with  (3.48)  and  (3.51),  in  (3.46),  to  obtain  an  equation  from  which  the  coeffi¬ 
cients  C2n,  n  >  1,  can  be  determined  by  comparing  like  powers  of  p.  Again  using  the  Mathematica  computer 
algebra  system,  we  find  the  coefficients  c-m ,  n  >  1,  to  be  given  in  terms  of  bo  by 


„  _  72n  ,2 
C2n  —  ySn  60. 


Tl  >  1, 


(3.58) 


where  the  72 „  are  purely  numerical  coefficients  given,  for  1  <  n  <  9,  by 


,  1  5  55  7 

72  =  1,  74  =  76  =  <j.  7s  =  j2’  710  =  6’ 


205 


17051 


712  108’  714  5292  ! 


2864485 
716  "  508032  ’ 


718  = 


103863265 
10287648  ' 


(3.59) 


The  remaining  constant  Co  is  determined  by  the  clamped  edge  boundary  condition  w  ( 1 )  =  0,  hence  from 
(3.56): 

OO  OO 

72 n 

b?? 


Co 


--E<*>  =  tfEiifc- 


(3.60) 


71—1 


Ei  uo 


71=1 


After  determining  Co,  the  dimensionless  axial  displacement  component  w(p)  is  given  by  the  power  se¬ 
ries  (3.56). 

Power  series  solutions  for  the  displacement  components  (as  well  as  the  stress  and  strain  components) 
having  the  proper  physical  units  can  be  obtained  by  restoring  the  scale  factors  introduced  in  equations 
(3.32),  (3.34),  and  (3.45).  For  example, 


u(R)  =  au(R/a )  =  \q2/3R 


-r  (1  —  va  )  +  bo  (1  —  va  )  +  bo  (2n  +  1  —  va  )  52n(-R/fl)2n 


71=1 


and 


w{ 


j(R)  =  aw(R/a)  =  aq1^3  w(R/a)  =  aq1/3  ^  C2n(R/a)2n 9 
and  the  apex  deflection  w0  =  iu(0)  is,  from  (3.62): 


71=0 


w0  =  ag1/3  c0. 


,  (3.61) 


(3.62) 


(3.63) 


3.3  Applications  to  Bulge  Testing 

The  last  equation  (3.63)  is  useful  for  the  analysis  of  data  obtained  from  bulge  testing  (see,  for  example,  [8, 
9,  10]  for  descriptions  of  such  tests),  a  method  under  consideration  for  experimental  determination  of  the 
intrinsic  coating  stress.  A  brief  discussion  of  such  an  analysis  will  follow,  but  first  we  compare  graphs  of 
apex  deflection  versus  pressure,  using  (3.63),  to  those  obtained  from  FE  analysis.  Figure  3  shows  the  results 
when  M  =  0,  i.e.,  for  a  zero  net  residual  stress  resultant.  The  coated  membrane  radius  here  is  a  =  3.0 
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Apex  Displacement  -vs-  Pressure 


Figure  3:  Comparisons  of  w0  versus  p,  using  equation  (3.63),  to  results  using  the  finite  element  method. 

These  plots  are  for  Af  =  0  (hence  also  r  =  r  =  0),  i.e.,  for  zero  net  residual  stress  in  the  coated  membrane 
laminate. 


inches  =  0.0762  m,  the  coating  and  membrane  thicknesses  are  hi  =  hc  =  1  pm  and  h2  =  hs  =  20  pm, 
respectively,  with  Ec  =  44.0  GPa  and  Es  =  2.2  GPa.  The  two  Poisson’s  ratios  have  been  assumed  the  same! 
i.e.,  vc  =  i 's  =  0.4.  FE  results  for  both  pinned  edge  and  clamped  edge  boundary  conditions  (see  §4.2.3  for 
the  definition  of  pinned  edge  boundary  conditions),  treating  the  coated  membrane  as  a  “plate”,  are  also 
shown  in  Figure  3.  It  is  clear  from  these  plots  that  the  results  of  the  theory  under  consideration  resemble 
more  closely  the  FE  results  for  a  plate  with  a  pinned  boundary  rather  than  one  with  a  clamped  boundary. 
We  interpret  this  to  mean  that  the  coated  membrane  behaves  more  like  a  plate  than  a  “true  membrane”,  i.e., 
the  bending  stiffness  of  the  coated  membrane  is  large  enough  in  this  example  that  a  geometrically  nonlinear 
plate  theory  is  probably  necessary  to  obtain  agreement  with  the  FE  clamped  edge  results. 

Returning  to  the  subject  of  bulging  tests,  we  recall  that  in  such  tests  the  apex  deflection  (“bulge”)  w0  is 
measured  as  the  pressure  load  p  is  varied.  The  analysis  of  data  obtained  in  this  way  often  relies  on  a  formula 
proposed  by  Beams  [8]  (see,  for  example,  [9,  p.  254],  cited  in  the  recent  book  by  Ohring  [10,  p.  716]).  Beams 
cites  for  this  formula  apparently  unpublished  work  by  Cabrera.  The  derivation  of  the  formula  is  not  given 
in  [8],  and  it  was  used  there  and  in  [9,  11]  only  to  analyze  thin  films  of  a  single  material,  not  a  laminate. 
A  similar  formula  was  given  in  [12]  for  applications  to  square  and  rectangular  films  of  a  single  material. 
These  formulas  seem  to  be  intended  for  applications  to  “true  membranes” ,  since  they  do  not  include  bending 
stiffness.  The  formula  given  by  Beams  [8],  expressed  in  our  notation  (Beams’  application  was  to  a  sinqle 
thin  film  of  gold),  is: 


__  ^  h  wp 
a  a 


To  + 


2  E 

3  1  -  v  a2)  ’ 


equation  (2)  of  [8], 


(3.64) 
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where  T0  is  referred  to  as  the  “tension”  by  Beams,  but  has  units  of  stress  and  corresponds  to  our  Af/h.  Since 
Poisson’s  ratio  for  gold  is  approximately  0.4,  for  Beams’  application  equation  (3.64)  can  be  rewritten  as 


p  =  4^A r  +  4.44£/i^f, 
a1  o4 


A f  =  hT0. 


(3.65) 


For  a  single  material,  in  which  case  E  =  E,  the  two  terms  of  the  Beams/Cabrera  formula  follow  as  an 
approximation  of  our  theory  by  investigating  the  apex  deflection  as  a  function  of  pressure  in  the  two  limiting 
cases  of  zero  residual  stress,  and  large  residual  stress.  For  zero  residual  stress  (Af  =  0),  equation  (3.63) 
applies,  where  cq  is  a  function  of  Poisson’s  ratio  only,  which  we  find  to  have  the  value  0.626  for  v  =  0.4,  i.e., 


Wo 


=  0.626  a  (|^)  /  )  AT  =  0,  v  =  0.4. 


(3.66) 


On  the  other  hand,  for  large  values  of  Af,  i.e.,  of  r,  the  coefficient  bo  in  (3.60)  is  found  to  be  large,  so  that 
co  «  l/&o-  In  this  case  it  also  follows  from  (3.55)  that  b0  «  r,  hence  Co  «  1/r  and  from  (3.63) 

q1/ 3 

wo  s=»  a-^-. 


But  from  (3.45),  r  =  4 r/g2/3,  which  yields 


Wofa4^q=4AfP’  ^ large’ 


(3.67) 


where  we  used  (3.45)  again  to  replace  q  and  r.  Thus,  for  large  A f,  w0  is  linear  with  pressure  p.  Equations 
(3.67)  and  (3.66)  yield  two  expressions  for  the  pressure  as  a  function  of  wo'- 

w0 


Pi 


:  4  — —  Af,  AHarge, 


and 


p2  '=  4.08  Eh 


> 


Af  =  0,  v  =  0.4, 


(3.68) 


(3.69) 


respectively.  The  Beams/Cabrera  formula  (3.65)  is  obtained  as  a  simple  linear  combination  of  these  expres¬ 
sions  for  pi  and  P2,  viz., 

(3.70) 


4.44 

”  =  +  408 


Since  this  formula  is  an  approximation,  it  is  expected  not  to  be  as  accurate  in  determining  the  intrinsic  stress 
as  one  determined  from  the  exact  theory.  A  more  systematic  derivation  of  Beams/ Cabrera-like  formulae  using 
perturbation  techniques  will  be  given  later  in  §5.2. 

Our  algorithm  for  using  the  theory  to  compute  the  coating  stress  Sc  in  a  single-layer  coating  from  bulge 
test  data  is  the  following.  For  a  given  pressure  p,  a  measurement  of  the  apex  displacement  w0  is  made.  With 
p  we  can  compute  q  from  (3.33),  and  then  from  (3.63)  we  have  the  corresponding  value  of  c0: 

CO  =  -3W.  (3.71) 


aq1/3 ' 

For  this  value  of  Co,  we  now  determine  bo  by  solving,  using  some  numerical  method,  equation  (3.60),  i.e. 

oo 

72  n 


n=l  °° 


(3.72) 


rather  than  (3.55).  This,  in  effect,  gives  us  bo  as  a  function  of  Co,  i.e.,  bo(co).  Substitution  of  this  value  of  bo 
in  (3.55)  then  gives  t  as  a  function  of  i>a  and  bo  (or,  equivalently,  Co): 


t(co,va)  =  b0 


1  + 


^2  (2n  +  1  -  VA  )  &2n 


(3.73) 
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X  -vs- 


c 


0 


Figure  4:  r  versus  cq  for  five  values  of  va- 


We  assume  that  the  geometrical  and  material  properties  of  the  coating  and  membrane  are  known,  so  that 
vA  can  be  calculated  from  its  definition  (3.33),  and  then  (3.73)  allows  us  to  compute  t(cd)  for  this  fixed 
value  of  va-  As  the  final  step,  the  computed  value  of  r  is  used  in  definitions  (3.45)  and  (3.33)  to  obtain  M, 
from  which  Sc  follows  according  to  definition  (1.1),  assuming  that  the  membrane  substrate  residual  stress 
Ss  is  known.  Figure  4  shows  the  graphs  of  f(co)  versus  Co  for  five  different  values  of  vA  (note  that  the  curve 
for  va  —  0.3,  restricted  to  the  interval  0.1  <  Co  ^  0.653,  is  the  inverse  of  the  graph  shown  in  Figure  6  of 
Campbell  [6]). 

In  closing  this  Section,  we  note  that  in  order  to  include  the  bending  stiffness  of  the  material,  Bonnotte, 
et  al  [13]  appended  an  additional  term  to  a  Beams/Cabrera-like  formula.  They  then  used  FE  analysis  to 
determine  two  fitting  constants  appearing  in  their  formula.  This  model  was  applied  later  by  Bahr,  et  al  [14] 
to  a  square,  five-layer  laminate  material,  by  simply  replacing  E  and  v  in  the  Bonnotte,  et  al  formula  with 
thickness-weighted  “composite”  values,  i.e., 

5  5  5 

E  =  E  hiEi/h>  v  =  hiVi/h>  h  =  Si  ki '  (3.74) 

«=1  «=1  t=l 

These  composite  values  appear  in  their  equation  (1)  only  in  the  combination  E/(  1  -  i/2),  which  they  refer 
to  as  the  composite  “biaxial  modulus”,  denoted  by  YB  in  Bonnotte,  et  al  [13],  and  by  Q  in  our  definition 
(1.8).  It  should  be  noted  in  this  regard  that  in  other  recent  literature,  e.g.,  [15,  16,  17],  the  term  “biaxial 
modulus”  is  used  instead  for  the  quantity  E/(  1  -  v)  (a  convention  which  we  shall  adopt  later,  denoting  it 
by  B).  For  a  Poisson’s  ratio  of  0.4  the  first  definition  gives  a  value  of  1.19  E,  while  the  second  gives  a  value 
of  1.67  E,  roughly  409c  higher  than  the  first.  In  Table  1  we  use  the  Bahr,  et  al  [14]  data  to  compare  values 
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Table  1:  Comparison  of  thickness- weighted  averages  with  values  occurring  in  the  geometrically  nonlinear 
membrane  laminate  theory. 


E  (GPa) 

V 

E/(l-u2)  (GPa) 

Thickness-weighted  averages 

145.6 

0.280 

158.0 

Membrane  laminate  theory 

145.8 

0.283 

158.5 

computed  using  the  thickness-weighted  averages  in  (3.74)  with  those  computed  using  the  definitions  (3.33) 
occurring  in  the  geometrically  nonlinear  membrane  laminate  theory.  For  all  practical  purposes,  the  values 
obtained  by  the  two  different  methods  axe  the  same  in  this  particular  example. 

Results  obtained  with  the  formulas  used  in  [13]  and  [14]  are  probably  not  comparable  to  results  using  our 
theory,  as  their  formulas  were  determined  empirically  by  fitting  to  data  obtained  on  rectangular  specimens, 
rather  than  circular  ones.  However,  ignoring  this  important  distinction  for  the  moment,  and  using  the  data 
given  in  [14],  we  compare  in  Figure  5  the  results  of  the  theory  with  the  results  shown  in  Figure  1  of  [14]. 
Note  that  their  data  is  for  a  multilayer  consisting  of  four  coating  material  layers  placed  on  a  (fifth)  silicon 

Pressure  -vs-  Apex  Displacement 
Sj  =  c.  =  1 17  MPa  in  each  coating  layer  i,  Diameter  =  1.45  mm 


Figure  5:  Comparison  of  membrane  laminate  theory,  equation  (3.63),  with  curve  generated  by  equation  (1) 
of  Bahr,  et  al  [14],  assuming  Si  =  ar  =  117  MPa  in  each  coating  layer  i. 


substrate  layer,  so  that  the  constants  of  our  theory,  defined  in  equations  (1.1)— (1.7),  involve  N  =  5  layers. 
For  the  theory,  we  took  the  radius  to  be  a  =  0.725  mm,  half  that  of  the  side  length  1.45  mm  of  their  square 
specimens  (in  [14]  this  side  length  is  also  denoted  by  a,  which  should  not  be  confused  with  our  radius  a). 
In  order  to  match  the  final  point  (3  x  10-5  m,1.4  x  105  Pa)  of  their  Figure  1  using  their  equation  (1),  we 
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-fn  ^d6  t^f*r  residuaI.  stress  °t  to  be  117  MPa,  although  they  mention  the  range  to  be  from  90  MPa 
to  110  MPa.  Our  theoretical  curve  was  generated  using  equation  (3.63),  assuming  that  S,  =  117  MPa  for 
each  layer  i  of  the  four-layer  coating,  and  S5  =  0  in  the  substrate.  The  resulting  comparison  is  reminiscent 
of  our  comparison  of  the  theory  and  FE  analysis  in  Figure  3,  i.e.,  for  a  given  pressure  the  theory  generally 
overestimates  the  corresponding  apex  deflection. 


4  Geometrically  Linear  Shell  Laminate  Problems 

For  shell  laminates,  the  leading  order  displacement  components  are  given  by  the  Kirchhoff-Love  expressions 


Ur  =  u  -  Zw,n,  Ue  =  v  -  Z  Uz  =  w, 

it 


(4.1) 


where  u  v,  and  w  are  functions  of  R  and  0  only.  For  a  geometrically  linear  shell  laminate,  the  strain  com¬ 
ponents  depend  linearly  on  these  components  and  their  partial  derivatives,  according  to  equations  (10  23  I)- 
(10.28.1)  of  Reference  [1]: 


e*e  =  j 


'  R  R 


=  e°RG  ~ 


£rr  uyR  +  r tRWtn  —  Z  kRR, 

-f  u 


«*>  =  ~ R -  -  Z(lT  +  ^B?)  E  -  Zke°- 

where  the  Z-independent  strains  and  curvatures  are  given  by 

=  *  f.,.  ,  "e-»  ,  r ,RWte\  ,  w?> 

Re  ~  2  V  1  +— R~  +  “r— J  ’  kRe  =  ~lt2  + 


w,fle 
R  ' 


€RR  =  U,R  +  r  RW>R,  kRR  =  W'RR, 


,0 

eee 


_  a,e  +  u  _w}R  wtee 

~  R  ’  ke6  =  -R-¥~W" 


(4.2) 

(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 


The  stress  resultants  and  couples  for  a  shell  laminate  are  given  by  equations  (10.39.1)— (10  44  I)  of  Refer¬ 
ence  [1]:  '  ’ 


Nr  =  Af  +  Ae°RR  +  A,e©0  +  BHrr  +  Bvk&Q,  (4 .8) 

Nrq  =  Ae  (°Re  +  Bq  kRQ ,  (4.9) 

Nq  =  Af  +  Ave°RR  +  yle@©  +  B„kRR  +  Bkee,  (4-10) 

Mr  =  -M  -  BeRR  -  But%e  -  DkRR  -  D„k&e,  (4.11) 

Mrq  =  -Bq  e°Re  -  De  kRe ,  (4.12) 

Me  =  -M  -  BveQRR  -  Bcqq  -  DvkRR  -  Dkee ,  (413) 
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but  note  that  the  Z-independent  strains  are  now  replaced  by  (4.5)-(4.7)  of  the  present  Section.  The  equi¬ 
librium  equations  for  a  shell  laminate  are  equations  (10.36 .1)-(10.38 .1)  of  Reference  [1],  which  are  similar  to 
those  for  a  membrane,  but  contain  additional  terms  involving  the  stress  couples: 

(■ RNr)iR  -  Ne  +  NR@,e  =  0,  (4.14) 

(R2Nrq)r  +  RNe,e  =  0,  (4.15) 

and 

jirT^TVfl  4-  ( RMr)  r  -  M@  +  r 

+  r ,rNr@  +  jg  (tfM*©)^  +  |  Me,©  +  (p  +  7o<?)  R  =  0.  (4.16) 

4.1  Reduction  to  an  Axisymmetric  System 

We  specialize  now  to  the  case  of  an  axisymmetric  system  by  assuming  that  none  of  the  variables  depend  on 
the  angular  coordinate  ©.  Thus,  all  terms  involving  partial  derivatives  with  respect  to  0  vanish,  leaving  the 
following  simplified  set  of  strain-displacement  relations: 


4© 

-  !(■•«- !)■ 

35- 

5o 

II 

o 

(4.17) 

eRR  = 

1  u,r  +  r,Rm,R , 

kRR  =  W,RR, 

(4.18) 

0  «  , 

€©©  —  ^  )  K©© 

W,R 

R  ’ 

(4.19) 

and  equilibrium  equations: 

(RNr)<r  -  Nq 

=  o, 

(4.20) 

{R2Nr@)r  = 

0, 

(4.21) 

[RT,rNr 

+  (RMr)'R  -  Me] 

+  (p  +  Jo9)R  =  o. 

,R 

(4.22) 

As  in  §3,  cf.  equations  (3.16)  and  (3.17),  equations  (4.21)  and  (4.22)  can  be  integrated  immediately.  And, 
in  order  for  the  results  to  be  regular  at  R  =  0,  we  must  again  set  the  integration  constants  to  zero,  which 
yields  the  results 


Nrq  =  0,  for  all  R.  (4.23) 

RT,r  Nr  +  (RMr)  r  -  Me  +  (p  +  7o  g)  ^  =  °-  (4.24) 

Since  Icrg  =  0,  we  have  from  (4.9)  and  (4.23)  that  e°Re  =  0  (since  Aq  is,  in  general,  nonzero).  Thus,  as  in 
§3,  it  follows  that  v(R)  =  no  R,  where  no  is  an  arbitrary  integration  constant;  if  v(Rq)  =  0  for  any  non-zero 
value  Ro  of  R,  then  v(R)  =  0  for  all  R.  Since  both  e°Re  =  0  and  kn@  —  0,  it  follows  from  (4.12)  that 

Mr©  =  0,  for  all  R.  (4-25) 
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The  remaining  system  of  equations  to  be  solved  is  then 


R^r,r  +  Nr  -  TV©  =  0, 


(4.26) 


where 


and 


T>  2 

R  MR'R  +  Mr  -  Me  +  R T,r  Nr  +  (p  +  7o  g)  —  =0, 


Nr  =  M  +  A  eRR  +  Av  €qq  +  B  kRR  +  Bv  kee , 
Ne  =  N  +  At,  e°RR  +  A  cqq  +  Bv  kRR  +  B  kee , 
Mr  =  -M  -  Be°RR  -  Bv€%e  -  DkRR  -  D„kee, 
Me  =  -M  -  Bve°RR  -  Bcqq  -  DvkRR  -  Dkee, 


e°RR  -  U,R  +  F,rw,r,  kRR  =  W'RR, 


(4.27) 

(4.28) 

(4.29) 

(4.30) 

(4.31) 

(4.32) 


c0  _  u  ,  _  W,R  , 

60©  5  ^©©  —  •  (4.33) 

Next,  as  in  §3.2,  we  transform  these  equations  to  dimensionless  forms  by  introducing  a  dimensionless 
coordinate  />,  dimensionless  displacement  components  u  and  w ,  and  a  dimensionless  initially  curved  reference 
surface  F,  defined  by 

R  ^  u  ^  w  ^  r 

p  =  t*  =  — ,  w  =  r  = 


a  a 

together  with  the  following  dimensionless  constants: 

^A  =  ^,  E  =  J-  (l  -  v%)  ,  T  =  ~~ 

A  h  K  A>’  Eh  '  Eha’ 

and  four  dimensionless  dependent  variables  xr,  x$,  yr,  and  ye  defined  by 


a 


M 

P  =  =- 


_  _  Nr- AT 
Eh 


_  Ne-M 
X°  “  Eh  ’ 


_  Mr  4*  M 


=  (P  +  7o  9)a 
q  Eh 


_  M©  +  M 


Eha 


ye  =  — = 


Eha 


We  also  introduce  dimensionless  curvatures  kt  and  kg  defined  by 

Kr  =  akRR  —  w 


w 


\pp 


kq  =  a  /c©©  =  — — , 
P 


and  four  new  dimensionless  constants 


b  = 


B 

Eha ’ 


Eha 


d  = 


£> 


£/ia2’ 


dy  —  — 


£/ia2’ 


(4.34) 


(4.35) 


(4.36) 


(4.37) 


(4.38) 


noting  that  the  strain  components  are  already  dimensionless,  and  have  the  following  forms  in  terms  of  the 
dimensionless  variables  defined  above: 


eRR  —  u,p  +  r<pWiP,  £00  —  — . 

P 


(4.39) 
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(4.40) 


The  strain  components  must  satisfy  a  modified  form  of  the  compatibility  condition  (3.40),  viz., 

Peee,p  “  eRR  “  e0©  ^,pwiP> 
and  there  is  an  additional  compatibility  condition  involving  the  dimensionless  curvatures: 

P  =  R't  K'd't  (4.41) 

which  is  easily  derived  from  the  definitions  (4.37). 

Substituting  these  dimensionless  quantities  into  the  equilibrium  equations  (4.26)  and  (4.27),  we  reduce 
them  to 


and 


#r,p  (  *Er  &6  )  — 

yr,p  +  -  (j/r  -  Ve )  +  f,p  (r  +  xr)  +  q  7.  =  °> 
P  1 


respectively.  The  constitutive  relations  (4.28)-(4.31)  take  the  forms 

Xr  =  1 - 2  (e°RR  +  ^e©©)  + 

A  “  V  a 


(4.42) 

(4.43) 

(4.44) 


xg  = 


(^A  ^rr  "l"  £©©)  4“  bv  Kr  +  bug, 


yr  —  ~bt°RR  bv  c@@  d  Kr  duKg, 
yo  — -  by  ^rR  6e©e  dv  Kr  dug. 

Equations  (4.44)  and  (4.45)  are  easily  solved  for  eRR  and  Cqq  in  terms  of  the  z’s  and  k’s: 

£ RR  —  Xr —  VAXg  P  Kr  pv  Kg, 

€00  —  Xg  V A  Xr  Pv  Kr  P  Kg  •> 


(4.45) 

(4.46) 

(4.47) 

(4.48) 

(4.49) 


where  we  have  introduced  two  new  constants: 


P  =  b  -  vAbv,  =  b„  -  vAb.  (4.50) 

We  note  here  that  the  dimensionless  radial  displacement  component  can  be  expressed  in  terms  of  xr  and  xg, 
and  Kr  and  Kg  (which  themselves  depend  only  on  derivatives  of  w),  by  combining  (4.39)  and  (4.49)  to  obtain 

u(p)  =  p  (xg  -  VA  xr  -  Kr  -  p  Kg)  .  (4-51) 

Tahing  the  derivative  with  respect  to  p  of  (4.49)  and  substituting  the  result,  together  with  (4.48)  and 
(4.49),  in  the  e-compatibility  relation  (4.40),  we  obtain 

p  (xg,p  -  UaXt,p  -  P„KT,p  -  P  Kg>p)  =  (1  +  Va)  (xr  -  Xg)  -  (P  -  P„)  {kt  ~  Kg)  -  T>pW>p. 

In  this  expression,  we  use  (4.41)  to  replace  (kt  -  Kg),  and  (4.42)  to  replace  ( xr  —  xg),  which  yields  the 
following  result  for  the  compatibility  condition: 

p  [ XT  +  Xg  -  pv  (Kr  +  Kg)\  p  +  ftPWlP  =  0.  (4.52) 
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Next,  we  substitute  (4.48)  and  (4.49)  in  (4.46)  and  (4.47)  to  obtain 

yr  =  ~Pxr  -  pvxe  -  5kt  -  6„i tg,  (4.53) 

y$  Pv  %T  P  Xg  Sy  Kr  8  K$ ,  (4.54) 

where  we  introduce  another  pair  of  constants: 

6  =  d  -  bf)  -  bvpv,  8„  =  dv  -  bpv  -  bvp.  (4.55) 

The  derivative  of  yT  with  respect  to  p  yields  from  (4.53): 

2/r,p  =  ~0%r,p  ~  Pt/Xgtp  <5  «r,p  —  8vKgtp.  (4.56) 

From  equations  (4.53),  (4.54),  and  (4.56)  we  obtain 


Vr,p  +  ~{yr  -  ye)  =  -PxT<p  -  Pvxe,p  -  SKrtP  -  8„  KgtP 

fr)(Xr  -  *e)  +  (6  -  Sv)(Kr  -  K0)}. 

In  this  expression,  we  again  use  (4.41)  to  replace  (/cr  -  Kg),  and  (4.42)  to  replace  {xr  -  x0),  which  simplifies 
it  to 

Vr,p  +  -  (yr  -  ye)  =  ~Pu  ( xr  +  xg)  p  -  6  (Kr  +  Kg)p.  (4.57) 

Substitution  of  (4.57)  into  (4.43)  yields  the  following  form  of  the  axial  equilibrium  equation: 


-ft,  (xr  +  Xg)  p  -  6  (nr  +  Kg  )  p  +  fp  (r  +  xr)  +  q  |  =  0.  (4.58) 

Equation  (4.42)  can  be  used  to  eliminate  xg  in  (4.52)  and  (4.58),  and  the  problem  we  are  left  with  is  to 
solve  the  coupled  differential  equations  (4.52)  and  (4.58)  for  w  and  xT.  The  radial  displacement  u  is  then 
determined  from  (4.51). 

Equations  (4.52)  and  (4.58)  differ  in  form  from  those  for  a  shell  of  a  single  material  only  by  the  terms 
involving  the  coefficient  pv.  From  its  definition,  given  in  (4.50),  we  note  that 


=  bu  -  vAb  =  =—  (uB  -  v A), 
hi  ha 

i.e.,  f}„  is  proportional  to  the  difference  in  extensional  and  coupling  Poisson’s  ratios.  This 
expanded,  yielding 


(4.59) 

expression  can  be 


Pv  = 


1 


AEha 


(. AB„  -  A„B) 


1  j  N  TV— 1  N 

~  AF.hr,  2  EE  E  hi  hjhk  [ Qi  {QjVj  Qkvk  )  —  Qi^i  ( Qj  —  Qk)]i 


AEha  2 


i=l  j=  1  k  =  j+ 1 


which  can  be  manipulated  to 


j  ^  Ar  — 1  TV— 1  TV 

A  E  h  a  i  2  EE  E  hi  hj  hk  Qi  Qk  ( Vi  Vfc ) 

i=i  j  =  i  t  =  j+i 

hNQN  ^  A  , 

+  — 2 —  E  E  hihk  \Qi  (vi  “  v*)  ~  Qk  (uk  -  vn)]  }• 

i  =  i  k  =  j+ 1 


(4.60) 
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This  constant  is  seen  to  be  a  sum  of  terms,  each  of  which  depends  upon  differences  in  the  Poisson’s  ratios  of 
the  various  material  layers.  For  a  single-layer  coating  of  the  membrane  (hence  N  =  2  layers  total),  we  have 
simply 

=  hchsQcQs.  {j/c_Ush  (4.61) 

2  AEa 

where  i  —  j  =  1  =  c  in  the  coating,  and  k  =  2  =  N  =  s  in  the  membrane  substrate.  The  total  thickness 
h  =  hc  +  hs  disappears  from  the  expression,  as  it  is  a  common  factor  in  both  numerator  and  denominator. 
Thus,  if  the  two  materials  have  the  same  Poisson’s  ratio,  /?„  vanishes  and  the  governing  equations  simplify 
considerably  (this  simplifying  assumption  was  made,  for  example,  by  Wittrick  [18]  in  his  work  on  the  stability 
of  bimetallic  thermostats). 


4.2  General  Solution  for  an  Initially  Parabolic  Coated  Membrane  Laminate 

Here,  we  consider  an  initially  parabolic  coated  membrane,  defined  by  equation  (1.19),  viz., 


r(-R)  =  §  (a2  -  R 2) , 

The  dimensionless  form  of  this  equation  can  be  written  as 

tip)  =  y  (i  -  p2)  - 

where  is  a  dimensionless  parameter  defined  by 


k0  =  ok  = 


2/' 


From  (4.63)  we  obtain 

F,p  —  —k0  Pi 

which  can  be  substituted  in  the  fundamental  equations  (4.52)  and  (4.58)  to  write  them  as 

[xr  +  Xe  -  Pv  («r  +  Ko)],p  ” 

-Pv  [xr  +  xe)  -  S  (/Cr  +  Kg  )  -  K0  p  (r  +  xT)  +  q  5  =  0. 


(4.62) 

(4.63) 

(4.64) 

(4.65) 

(4.66) 

(4.67) 


This  pair  of  equations  can  be  uncoupled  to  obtain  a  single  differential  equation  for  w  as  follows.  Substituting 
for  (xr  +  xs)  p  from  (4.66)  into  (4.67)  yields 


-A  (/cr  +  Ke)  p  -  PvK0w,p  -  K0p  (r  +  xr)  +  q  -  =  0, 
where  we  have  introduced  a  new  constant  A  defined  by 

A  =  S  +  Pl  =  d  -  (1  -  v\)  b\ 

Solving  (4.68)  for  xr  yields  (for  ^  0) 


Xr  =  — T  — 


A(*r  +  **)* 


w.o  q 

P  2 


from  which 


xr,p  — 


(Kr  +  Ke  )tPP  (Kr  +  He  )„ 


a  (Ops.  _ 
v  V  p  p1  )' 


(4.68) 

(4.69) 

(4.70) 

(4.71) 
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FYom  equations  (4.42),  (4.70),  and  (4.71)  we  have 


A 

xe  —  Pxr,p  +  xr  =  “7"  —  —  («r  -f  Kq  )  pp  ~  + 

K/\  ’ 


Q 

2k0' 


hence  (4.70)  and  (4.72)  yield 


xr  +  xg  =  -2r  +  — —  — 

Kq  Kq 


We  next  observe  that  equation  (4.66)  can  be  integrated  immediately  to  obtain 

xr  +  xg  =  /?„  (kt  +  Ke)  +  k0w  +  C3, 

where  C3  is  an  arbitrary  integration  constant.  Substituting  from  (4.73)  into  (4.74)  yields 


(4.72) 


(4.73) 


(4.74) 


A 

Ka 


(«r  +  *6  )pp  + 


(« r  +  Kfi  ) 


+  Pi 


tc  p 

pp  H - —  +  Kr  +  Ke 

P 


j  +  K0 


w  +  2r - h  C3  =  0.  (4.75) 


V2w  =  uJ.pp  + 


Recalling  the  definition  of  the  p-dependent  part  of  the  Laplacian  operator  in  cylindrical  coordinates,  viz., 

w,p 

-f’  (4.76) 

(4.77) 

(4.78) 


we  have 


and  also 


w  1 


+  K0  =  +  —  =  V2u5, 

p 


V2  (V2t2)  =  V4ui  =  (V2u>)  pp  + 


These  observations  allow  us  to  write  equation  (4.75)  more  compactly  as 

+  2 W  +  4™  =  X  («  -  2tk0  -  Csk0). 


(4.79) 


The  complete  solution  of  the  linear  differential  equation  (4.79)  has  the  form  w(R)  =  wh(R)  +  wp(R ),  where 
wh  (R)  is  the  general  solution  of  the  homogeneous  equation 


V*wh  +  2  + 


(4.80) 


and  ihp(R)  is  a  particular  solution  of  (4.79).  A  suitable  particular  solution  in  this  case  is  the  constant  function 


Wp(R)  =  ~2  (Q  ~  2r«0  -  C3K0), 


(4.81) 


which  involves  the  as  yet  unknown  integration  constant  C3.  The  homogeneous  fourth-order  linear  differential 
equation  (4.80)  has  at  most  four  linearly  independent  solutions.  We  note  that  if  ip(p)  is  an  eigenfunction  of 
the  p-dependent  Laplacian  operator  defined  by  (4.76),  corresponding  to  an  eigenvalue  -m2,  i.e.,  if 


V2V>  =  -m2ip,  (4.82) 

then  VV  =  mV,  hence  tp  is  also  a  solution  of  (4.80)  ifm  satisfies  the  fourth-degree  polynomial  equation 

m4  +  2/cfm2  +  k\  =  0,  (4.83) 
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where  we  have  introduced  two  constants  &2  and  ki  defined  by 


kn 


(4.84) 


and  remark  that  j}v  can  be  either  positive  or  negative,  while  A  is  always  positive. 
(4.83)  are  easily  found  to  be  , _ 


m 


The  four  solutions  of 


(4.85) 


The  forms  of  these  solutions  to  be  used  depends  on  the  relation  between  k\  and  k\ .  We  find  for  typical 
parameter  values  of  the  material  constants  and  the  geometry  that  kf  >  k\,  so  that  (4.85)  may  be  written  as 


m  =  ±)J~k%  ±iyj(k\  4-  kl){k\  -  k\)  . 

This  can  be  put  into  a  more  convenient  form  by  introducing  two  new  (real)  constants  defined  by 

k  =  \J  +  fcf)  ,  e  =  \J -  fc|)  i > 

which  can  be  inverted  to  yield  k\  and  k\  in  terms  of  e  and  k: 

k2  =  k2-e2,  k\  =  k  2+e2. 

The  complex  eigenvalues  (4.86)  take  the  following  simple  forms  in  terms  of  e  and  k : 

m  =  ±{e  ±  ik) . 


(4.86) 


(4.87) 


(4.88) 

(4.89) 


The  four  distinct  eigenvalues  mj,  j  =  1, 2,3,4,  are  thus  given  by 


mi  =  e  4-  ik , 
m2  =  e  -  ik  =  ml, 

771-3  —  —  (e  T  zfc)  =  -mi, 

7724  =  —  (e  —  ifc)  =  m^  =  —m2  =  —  m*, 


(4.90) 


where  an  asterisk  denotes  complex  conjugation.  These  correspond  to  eigenfunctions  'ipj  (p)  satisfying 

VVj  =  -ntfipj,  j  =  1, 2,3,4.  (4.91) 

Substituting  for  the  Laplacian  operator,  we  obtain  the  following  second-order  differential  equation  for  each  ipj : 

+  J  +  m2jiPj  =  0.  (4.92) 

Introducing  new  independent  (complex)  variables  =  rrij  p,  the  jth  second-order  equation  can  be  written 
as  , 

<Pj+T  +  =  °>  (4‘93) 

S  j 

where  a  prime  always  denotes  the  derivative  of  a  function  with  respect  to  its  argument,  in  this  case  and 
=  ^j(mj  p).  Equation  (4.93)  is  Bessel’s  equation  of  zero  order,  whose  solution  regular  at  the  origin 
is  the  zero-order  Bessel  function  J0(fj).  The  general  solution  of  our  homogeneous  equation  (4.80)  is  then  a 
linear  combination  of  the  four  eigenfunctions  Jo(£j): 

Wh(p)  —  Ci  Jo  (mi  p)  4-  C2  Jo(m2p)  4-  C3«7o(m3p)  4-  C4  Jo(m4p),  ,  ^ 

—  Ci  Jo(mip)  4-  C2  Jo(m*p)  4-  C3Jo(—mip)  -(-  C4  Jo(— m^p), 
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where  the  arbitrary  constants  Cj  are,  in  general,  complex.  Using  the  fact  that  J0  is  even,  i.e.,  J0(-z)  =  J0(z) 
this  homogeneous  solution  reduces  to  ’ 


Whip)  =  CiJo(mip)  +  C2Jo{mlp),  (4.95) 

where  Ci  =  C\  +  C3  and  C2  =  C2  4-  C4  are  new  arbitrary  complex  constants.  From  its  definition  in  (4.90), 
we  can  write  the  complex  number  mi  in  polar  form  as 


mi  =  e  +  ik  =  \/e2  +  fc2  e,e‘  =  ki  e* 6l , 


where  we  also  used  the  second  equation  of  (4.88).  Equation  (4.95)  is  thus  equivalent  to 

Mp)  =  CyJoie^'hp)  +  C2J0(e-i9'  klP). 

Now,  the  ratio  k/e  defining  the  angle  can  be  rewritten  using  (4.87)  and  (4.84)  as 


k  _  ll  +  /?„A-V2 

e  “  V  1  “  &A-1/2- 


(4.96) 


(4.97) 


(4.98) 


We  find  that  for  the  material  and  geometric  parameters  of  interest  here,  it  is  in  fact  true  that  |/?„  A-1/2|  is 
always  much  less  than  1  (in  a  typical  case,  this  product  is  on  the  order  of  1  x  10~8 ).  We  thus  approximate 


which  brings  (4.97)  to  the  form 


0,  ~  — 
"i  ~  4> 


(4.99) 


Whip)  —  Ci Joietn/4kip)  +  C2J0ie  t7!^kip).  (4.100) 

This  can  be  written  in  terms  of  Kelvin’s  her  and  bei  functions  (see  [19,  pp.  379-385],  for  details  of  the 
properties  of  the  Kelvin  functions)  as 


Whip)  —  Ci  [ber(fcip)  +  zbei(fcip)]  +  C2  [ber(fcip)  -  ibei(Arip)] ,  (4.101) 

Since  Kelvin’s  functions  are  real-valued,  and  wh(p)  must  be  real,  it  follows  that  the  complex  constants  must 
in  fact  be  complex  conjugates:  C2  =  C\.  Thus,  the  general  solution  for  w(p)  takes  the  form 


wip)  =  whip)  +  Wpip)  =  Ci ber(fci p)  +  C^beiihp)  +  \  (q  -  2tk0  -  C3k0),  (4.102) 

K0 


where  C\  _  C\  +  C\  and  C2  =  i(Ci  -  Cx)  are  now  arbitrary  real  constants  to  be  determined  from  the 
boundary  conditions. 

The  dimensionless  radial  displacement  u  defined  by  (4.51)  can  be  expressed  in  terms  of  w  and  its  deriva¬ 
tives  as  follows.  First,  we  use  (4.70)  and  (4.72)  to  replace  xr  and  xe  in  (4.51),  yielding 


From  (4.79)  and  (4.78)  we  obtain  the  following  expression  for  (V2u5)  ^ : 
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which  is  substituted  for  (V2w;)  pp  in  the  previous  equation  to  obtain  the  desired  result,  viz., 


w(p)  =  -  2^-^  (1  +  VA  )  +  C3  p  +  (1  +  Va)  — 

+  [Pv  (2  +  va)  ~  P]  w,p  +  K0pw.  (4.103) 

To  get  the  derivatives  of  w(p)  appearing  in  (4.103),  in  terms  of  Kelvin  functions,  we  note  that 

w,P(p)  =  fci  [C'iber'(fcip)  +  C2bei'(fcijo)]  ,  (4.104) 

wPP(p)  =  k\  [C1bev"(klP)  +  C2bei"(klP)] ,  (4.105) 

where  the  prime  on  the  Kelvin  functions  here  denotes  a  derivative  with  respect  to  the  function  argument  x  = 
kip.  Prom  [19],  Equations  (9.9.16), 


hence 


ber'(x)  =  —f=[  beri  (a:)  +  beii  (x)] , 

bei'(x)  =  -^=  [-beri(ar)  +  beii(x)] , 
V2 


ber"(a;)  =  -^=  [  beri  (a;)  +  beii  (*)] . 
bei"(x)  =  -^=  [-beri  (a;)  +  beii  (a;)]  , 


(4.106) 

(4.107) 

(4.108) 

(4.109) 


For  the  derivatives  of  the  Kelvin  functions  of  order  1,  we  have  from  Equations  (9.9.14)  of  [19]  the  following 
identities: 


ber'^x)  =  —  [ber2(x)  +  bei2(x)  -  ber(x)  -  bei(x)] , 


beij  (x) 


V2 


[bei2(x)  -  ber2(x)  -  bei(x)  +  ber(x)] , 


(4.110) 

(4.111) 


and 

r2 

ber2(x)  = - [beri(x)  —  beii(x)]  —  ber(x), 

x 

r2 

bei2(x)  = - [beii  (a)  +  beri  (a?)]  -  bei(x). 

x 

Applying  the  last  four  identities  to  (4.108)  and  (4.109)  yields 

ber"(x)  =  --ber'(x)  -  bei(x), 
x 

bei"(x)  =  --bei'(x)  +  ber(x). 
x 

Substituting  these  results  in  (4.105),  we  obtain 

w,Pp(p)  =  [Ciber^fcip)  +  C2bei'(fcip)]  4-  k\  [-C\  bei(fcip)  +  C2ber(fc1p)] 


w 


»p 


+  k\  [-Cl  bei(fcip)  +  C2  ber(A?ip)] 


(4.112) 

(4.113) 

(4.114) 

(4.115) 

(4.116) 

(4.117) 
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where  we  used  (4.104)  in  (4.116)  to  get  (4.117),  and  recall  that  x  =  klP.  Prom  (4.117)  follows  a  useful  form 
of  the  p-dependent  Laplacian  operator  acting  on  w: 


V2w  =  wiPP  +  ^  =  k\  [-Cibei(A:1p)  +  C2ber(A:ip)] .  (4.118) 

The  derivative  of  the  last  expression,  which  is  required  in  (4.103),  can  be  written  as 

(V2w)  p  =  kf  [-Ci  bei'(fcip)  +  C2  ber'^p)]  .  (4.119) 


Substituting  (4.118),  (4.119),  and  (4.102)  in  (4.103)  yields  the  general  solution  for  G(p)  in  terms  of  Kelvin 
functions  (the  terms  containing  C3  conveniently  add  to  zero): 


u{p)  =  Ulfcf  [-C1bei'(fc1p)  +  C'2ber'(fc1p)]  +  u2  h  [Ciber'^jp)  +  C2bei'(fcip)] 
+  Kop  [Ciber(fcip)  +  C2bei(A:1p)]  -  (1  -  vA)  (r  -  ^  p, 

where  we  have  introduced  new  constants  uY  and  u2  defined  by 


(4.120) 


,1  \  A 

ux  =  {  \  +  vA)  —  (4.121) 

Kq 

U2  =  pv  ( 2  +  vA )  -  p .  (4.122) 


Before  considering  specific  boundary  conditions,  we  note  that  for  any  such  conditions  that  include  a 
specification  w(p0)  =  0  at  some  point  p  =  p0  (where  p0  is  typically  either  0  or  1),  we  will  have  from  (4.102) 
a  condition  of  the  form: 


0  =  Ciber(fe!p0)  +  C'2bei(fc1p0)  +  \  ( q  -  2tkd  -  C3k0),  (4.123) 

K0 

■which  can  be  used  to  replace  the  constant  particular  solution  in  (4.102).  The  solution  for  uj(p)  is  thus  given, 
after  applying  a  boundary  condition  of  this  type,  by 

w{p)  =  Cl  [ber(fcip)  -  ber(fcipo)]  +  C2  [bei(fcjp)  -  bei(fcip0)].  (4.124) 


4.2.1  Computing  the  Kelvin  Functions 

For  arguments  of  the  Kelvin  functions  satisfying  kip  >  8  we  use  the  asymptotic  forms  of  the  Kelvin  functions 
and  their  first  derivatives,  derivable  from  material  given  in  Abramowitz  and  Stegun  [19]  (in  particular,  their 
Equations  (9.10.1)  and  (9.10.2)  on  p.  381,  together  with  Equations  (9.9.16)  on  p.  380),  viz., 


ber(kip)  ~  F(kip)  cos  ^!p/\/2  -  7r/8), 

(4.125) 

bei(fcip)  ~  F{kip)  sin  ^p/%/2  -  7t/8^ , 

(4.126) 

ber'^p)  ~  F(kxp)  cos  ^ip/\/2  +  7r/8)> 

(4.127) 

bei'(fcip)  ~  F(kip)  sin  (kip/y/2  +  7t/8), 

(4.128) 

noting  that  each  approximation  contains  a  common  multiplicative  factor 

V27r/ci  p 

(4.129) 
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The  sometimes  large  exponential  factor  F{k\p)  may  be  alternately  removed  and  inserted  in  constants  (defined 
later)  associated  with  various  boundary  value  problems,  in  such  a  way  that  it  will  ‘divide  out’,  to  avoid 
overflow  problems  in  the  computation. 

For  —8  <  kip  <  8,  we  use  the  polynomial  approximations  given  in  [19,  §9.11,  p.  384],  i.e., 

ber (kip)  =  1.0  -  64. Ox4  +  113.77778x8  -  32.36346x12 

+  2.64191a;16  -  0.08350x20  +  0.00123a;24  -  O.OOOOlx28,  (4.130) 


bei(fcxp)  =  16. Ox2  -  113.77778x6  +  72.81778a;10  -  10.56766a;14 

4-  0.52186a;18  -  0.01104a;22  +  0.00011a;26,  (4.131) 


her  (hp)  =  8a;  (-4.0a;2  +  14.22222a;6  -  6.06815a;10  +  0.66048a;14  -  0.02609a;18  +  0.00046a;22),  (4.132) 

bei'(fcip)  =  8a;  (0.5  —  10.66667a;4  +  11.37778a;8  —  2.31167a;12 

+  1.14677a;16  -  0.00379a;20  +  0.00005a;24),  (4.133) 

where  x  =  kip/S  on  the  right-hand  sides  of  these  formulas. 


4.2.2  Free  Edge,  Simply-Supported  at  the  Center 

The  first  boundary  value  problem  we  consider  is  that  of  an  initially  parabolic  coated  membrane  laminate 
with  a  free  edge,  requiring  that  Nr(cl)  =  MR(a)  =  0,  and  simply-supported  at  its  center,  i.e.,  w( 0)  =  0.  In 
terms  of  our  dimensionless  quantities  these  translate  to  xr(l)  =  —r ,  yr(  1)  =  p ,  and  w( 0)  =  0 ,  respectively. 
The  boundary  condition  w( 0)  =  0  yields  from  (4.124)  (with  p0  =  0  ): 

w(p)  =  Ci  [bev(kip)  -  1]  4-  C2bei(A;ip),  (4.134) 

since  ber(0)  =  1  and  bei(0)  =  0.  The  remaining  two  boundary  conditions,  xr  (1)  =  —r  and  yr(  1)  =  /x,  provide 
two  linear  algebraic  equations  to  be  solved  for  the  unknown  coefficients  C\  and  C2.  The  construction  of 
these  equations  involves  computing  the  strains  and  curvatures  defined  in  (4.39)  and  (4.37),  viz., 

U  7/} 

eRR  ~  «,/>  “  «0  PW,p,  €00  =  «r  =  W.pp,  ■  (4.135) 

where  we  replaced  fiP  =  -kq  p  in  the  first  equation  of  (4.135).  All  the  relevant  quantities  have  been 
computed,  except  for  utP.  From  (4.120)  we  obtain 


utP  =  U\  k\  j^-Cibei’^fcip)  +  C2ber  (fcip)]  +  u2fc2  [Ciber  (kip)  +  C^bei  (fcip)] 

+  K0pwtP  +  K0  [C'iber(fcip)  +  C2bei(fcip)]  -  (1  -  vA )  (r  -  ^)> 

or,  using  (4.114)  and  (4.115)  to  eliminate  the  second  derivatives  of  the  Kelvin  functions: 

utp  =  ui~  [C'ibei'(fci.p)  —  C2ber'(fcip)]  -  uifcf  [C’iber(fcip)  +  C2bei(fcip)] 

—  u2  —  [Ciber'^ip)  +  C^bei'tfcip)]  -  u2k\  [Cibei(fcip)  —  C2ber(fcip)]  (4.136) 

P 

+  K0pwtP  +  Kq  [C7iber(A:ip)  +  C2bei(kip)]  -  (1  -  vA)  (t  -  )  . 
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Using  these  results,  the  strains  and  curvatures  are  given  by 
k3 

t°RR  =  ui-j  [Cibei'fap)  -  Caber'fop)]  -  mfc4  [Ciberfop)  +  C2bei(Jfeip)] 

k 

-  U2  ~  [Ciber'^p)  +  C2bei'{k1p)}  -  u2k\  [Cjbeifop)  -  C2ber(^p)] 

+  K0  [Ciber(fcip)  +  ^bei^p)]  -  ( 1  -  vA )  ^  ,  (4.137) 

k3  Jr 

€ee  ~  Ul~^  [-Ci bei^Aqp)  +  Caber' (Aqp)]  +  u2 [C'iber'(fcip)  +  C2bei'(fcip)] 

+  k0  [Ciber(fcjp)  +  C2bei(*ip)]  -  ( 1  -  vA )  )  ,  (4.138) 


Kr  =  “7  [Ciber'(M  +  C2bei'(M]  +  k\  [-Cibei(fc1p)  +  C2ber(jk1p)] 
K°  ~  ~p  [Ciber^Aqp)  +  C2bei'(fcip)]  . 


(4.139) 

(4.140) 


These  expressions  must  be  evaluated  at  p  =  1,  and  then  substituted  in  the  constitutive  relations  (4.44)  and 
(4.46)  with  xT  =  -r  and  yr  =  p,  i.e., 

eRR  +  VA  C©e  +  (l  -  Va)  (bltr  +  bv  Ke)  +  (l  -  v\)  r  =  0.  (4.141) 

be°RR  +  6ree0  +  d/cr  +  dv  k8  +  p  =  0,  (4.142) 

The  resulting  system  of  two  equations  for  C\  and  C2  can  be  written  as 

snCj  +  si2C2  +  (1  -  i/j)  =  0,  (4.143) 

S21C1  +  S22C2  +  {0  +  0„)  2^-  +  p  -  (0  +  /?„)  t  =  0,  (4.144) 

whose  solutions  are  easily  found  to  be 

Cl  =  { [(1  -  ”a)  s22  ~  {P  +  0V)  S12]  2^ - S12  [p  -  (p  +  0„)  t]|,  (4.145) 

°2  =  |[(1“*'i).®2i  -  (P  +  P»)  «n]  -  Sn[p  -  {P  +  Pu)  T]|,  (4.146) 

where  |s|  =  «ns22  -  si2s2i  is  the  determinant  of  the  2  x  2  coefficient  matrix.  The  matrix  elements  are  rather 
complicated  expressions  involving  the  Kelvin  functions  and  material /geometrical  parameters.  They  can  be 
reduced  to  the  following  forms: 

«ii  =  [/co  (1  +  vA)  -  k\u  1]  ber(Aq)  -  k\  [u2  +  b  (l  -  v ’j|)]  bei(fci) 

+  k\ux  (1  -  vA)  bei'(fci)  +  ki(l  -  vA)  (P„  -  0  -  u2)  bet'fa),  (4.147) 

S12  =  [/c0  (1  +  vA)  -  /c4ui]  bei(fca)  +  kf  [u2  +  b  (l  -  v%)]  ber(fci) 

-  k\ux  (1  -  vA)  ber'(fci)  +  fci(l  -  vA)  {0V  -  0  -  u2)  bei'(fci),  (4.148) 
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521  =  [K0(b  +  bp)  -bkfui]  ber(fci)  -  h\{d  +  bu2)  bei(fci) 

+  kfui  ( b  -  bv)  bei'(fci)  -  fa  [d  -  dv  +  (b  -  bu)  u2]  ber'(fci),  (4.149) 

522  =  [kq  (b  +  bv)  -  b  k*ui]  bei(fci)  -1-  k\(d  4-  bu2)  ber(Ai) 

-  k\ux  (b  -  bu)  ber'(fci)  -  fa  [d  —  d„  +  (6  —  bv)  u2]  bei(fa).  (4.150) 

In  the  range  fa p  >  8  the  asymptotic  approximations  (4.125)-(4.128)  of  the  Kelvin  functions  contain  a 
common  exponential  factor  F(fap ),  having  the  value  F(fa)  at  the  edge  p  =  1.  Here,  we  see  that  in  the  same 
range  F(fa)  is  also  a  common  factor  of  each  matrix  element  hence  its  square  is  a  factor  of  the  determinant 
\s\.  It  then  follows  that  the  integration  constants  C\  and  C2  are  proportional  to  the  reciprocal  of  F(fa).  Since 
both  u  and  w  contain  products  of  Ci  and  C2  with  the  Kelvin  functions,  we  can  factor  from  each  occurrence  of 
these  functions  a  term  of  the  form  F(fap),  so  that  a  common  ratio  F(fap)/F(fa)  =  exp  [(fa/y/2  )(p  -  1)]/ yfp 
can  be  factored  from  these  terms.  Thus,  for  values  of  fap  >  8,  we  replace  the  Kelvin  functions  by  the 
trigonometric  parts  of  their  asymptotic  expansions,  and  factor  out  this  common  ratio  wherever  possible  to 
avoid  computational  problems  that  may  occur  for  large  values  of  the  exponential  function. 

It  is  rather  impractical  to  consider  a  pressure  difference  between  the  faces  of  a  coated  membrane  laminate 
satisfying  free  edge  boundary  conditions.  We  thus  set  the  pressure  difference  p  =  0,  hence  the  dimensionless 
constant  q  depends  only  on  the  gravitational  field.  Ignoring  the  effects  of  gravity,  we  set  q  =  0.  The  solutions 
(4.145)  and  (4.146)  for  C\  and  C2  are  then  observed  to  contain  the  common  factor 

p-  (/?  4-  ftp)  r,  (4.151) 

and  this  is  the  sole  dependence  of  C\  and  C2  on  the  intrinsic  stress  loads  Si  (occurring  in  p  and  r ).  We  can 
thus  write  the  solutions  in  (4.134)  and  (4.120)  for  w(p)  and  u(p)  (with  q  =  0)  as 

w(p)  =  [p-  (0  +  Mt]w(p),  (4.152) 

u(p)  =  [p  -  (l 3  +  ftp)  t)  u(p)  -  T  ( 1  -  I/A )  Pi  (4.153) 

where  w(p)  and  u(p)  are  functions  that  do  not  depend  on  the  intrinsic  stress  loads.  From  (4.152)  we  see 
that  the  axial  displacement  will  vanish  if  we  can  choose  the  intrinsic  stress  loads  and  geometrical/material 
parameters  to  satisfy  p  -  (ft  +  pv)  r  =  0 .  In  order  for  the  radial  displacement  to  vanish  as  well,  we  must 
also  choose  them  such  that  r  =  0.  This,  in  turn,  requires  p  =  0  for  the  first  condition  to  be  satisfied.  Thus, 
for  free  edge  boundary  conditions,  and  ignoring  the  effects  of  gravity,  the  necessary  and  sufficient  conditions 
for  there  to  be  no  displacement  from  the  initial  paraboloidal  shape  upon  removal  from  the  mold  are  simply 

r  =  0  and  p  =  0  =>  J\f  =  0  and  A4  =  0,  (4.154) 


that  is,  the  net  intrinsic  stress  resultant  and  net  intrinsic  stress  couple  must  both  be  zero.  Note  that  for  a 
single  coating  on  the  membrane,  the  condition  p  -  (ft  +  fiv)  r  =  0 ,  which  is  sufficient  to  insure  zero  axial 
displacement,  reduces  to 


Sc  -  B  Ss  —  0, 


where 


j2  —  Qc  (1  +  Vc) 

~  Q,-(l  +  i/s) 


(*)/(*)■ 


(4.155) 


and  we  made  use  of  the  first  definition  in  (1.8).  If  the  residual  stresses  are  thermally  induced,  so  that 


Si  =  - 


Ei 

1  ~n 


Oii  AT, 


where  a*  is  the  CTE  of  either  the  coating  (i  =  1  =  c)  or  membrane  substrate  (i  =  2  =  s),  then 

Sc  -  BSs  =  -  (3^7)  («c  -  a.) 


(4.156) 


(4.157) 
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which  vanishes  only  if  the  membrane  and  coating  have  the  same  CTE. 

Figures  6  and  7,  which  follow,  compare  our  solutions  for  the  axial  and  radial  displacements  with  geomet¬ 
rically  linear  finite  element  (FE)  solutions  of  the  same  problem.  The  model  considered  here  is  an  initially 
parabolic  membrane  with  a  single  coating.  It  has  a  radius  of  a  =  10  cm  =  0.1  m,  and  an  /-number  of  2, 
i.e.,  F&  —  2,  corresponding  to  k  =  1.25/m,  see  equation  (1.18).  The  coating  has  thickness  h\  =  hc  =  lpm, 
and  the  membrane  thickness  is  h2  =  h8  =  20pm.  The  coating  modulus  is  taken  to  be  Ec  =  44.0  GPa,  and 
the  membrane  modulus  Es  =  2.2  GPa.  We  assume  the  Poisson’s  ratios  of  the  coating  and  membrane  to  be 
the  same,  viz.,  vc~v8-  0.4.  The  assumed  intrinsic  stresses  in  this  example  are  small:  Sc  =  -5  KPa  and 
S8  =  20  Pa. 


u(R)  -vs-  R  Comparison 

Simply-Supported  Center  Boundary  Conditions  (Sc  =  -5  KPa,  S$  =  20  Pa) 


Figure  6:  Comparison  of  theory  to  geometrically  linear  FE  ^-displacement  results  (free  edge/simply- 
supported  at  the  center). 


4.2.3  Pinned,  or  Hinged,  Edge 

The  boundary  conditions  for  a  coated  membrane  with  pinned  (or  hinged)  edge  are  w(a)  =  0,  u(a)  =  0,  and 
Mr (g)  0.  In  terms  of  the  dimensionless  coordinate  p,  we  have  u;(l)  =  0,  so  that  po  =  1  in  equation 
(4.124),  2(1)  =  0  and  yr{  1)  =  /i.  The  first  condition  yields  from  (4.124): 

™(p)  =  Ci  [ber(fcip)  -  ber(fci)]  +  C2  [bei(fcip)  -  bei(Aq)].  (4.158) 

The  solution  for  u(p)  is  given  by  (4.120),  except  that  the  coefficients  Ci  and  C2  are  not  necessarily  the  same 
as  those  in  (4.120),  as  we  shall  soon  see. 

The  remaining  two  boundary  conditions,  2(1)  =  0  and  yr(  1)  =  p,  again  provide  a  system  of  two  equations 
for  the  constants  C\  and  C2.  We  note  that  since  2(p)  and  w,p{p)  have  precisely  the  same  forms  as  in  the 
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w(R)  -vs-  R  Comparison 

Simply -Supported  Center  Boundary  Conditions  (Sc  =  -5  KPa,  Sg  =  20  Pa) 


Figure  7:  Comparison  of  theory  to  geometrically  linear  FE  ^-displacement  results  (free  edge/simply- 
supported  at  the  center). 


simply-supported  center  problem,  yr(  1)  =  p  duplicates  equation  (4.142): 

b  ^jir  "b  by  “h  d  K,r  -f-  dv  kq  -f-  /i  =  0.  (4.159) 


The  equation  u(  1)  =  0,  however,  is  new  and  this,  together  with  (4.159),  leads  to  the  following  system  of 
equations: 


PiiCi  +  P12C2  -  (1  —  va)  “  0,  (4.160) 

P21C1  4-  P22C2  +  /i  —  (/?  +  Pv)  —  =  (4.161) 


where  P21  =  $21  and  p22  =  $22,  given  in  equations  (4.149)  and  (4.150).  The  matrix  elements  pn  and  pn 
have  the  comparatively  simpler  forms 


pn  =  Kober(ki)  —  kf  u±  bei'(ki)  H-  k\  ii2  ber'(A;i), 

P12  =  Kobei(fci)  +  ui  ber'(fci)  +  k\  u2  bei'(fci). 
Solving  (4.160)  and  (4.161)  for  C\  and  C2,  we  obtain 


Ci 

C2 


=  (1  -  OA)(r  -  J-)  -  Pn  [w  +  M  (r  -  £)  -  M  j , 

“'H{Bl(1“‘'J>(T_i)  ~Pn  [(',  +  w(T_i)  ■'•]}• 


(4.162) 

(4.163) 

(4.164) 

(4.165) 
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where  |p|  =  P11P22  -  P12P21  is  the  determinant  of  the  new  coefficient  matrix. 

Using  the  same  model  described  for  the  free  edge  problem,  Figures  8  and  9  compare  our  pinned  edge 
solutions  for  the  axial  and  radial  displacements  with  the  geometrically  linear  FE  solutions. 


u(R)  -vs-  R  Comparison 

Pinned  (Hinged)  Edge  Boundary  Conditions  (Sc  =  -5  KPa,  Ss  =  20  Pa) 


Figure  8:  Comparison  of  theory  to  geometrically  linear  FE  u-displacement  results  (pinned  edge). 


4.2.4  Rigidly  Clamped  Edge,  and  Coating  Stress  Prescriptions  for  Maintaining  an  Initially 
Parabolic  Shape 

A  near  net-shape  coated  membrane  used  as  the  primary  mirror  of  a  telescope  will  likely  be  attached  to  a 
rigid  circular  boundary,  hence  the  boundary  conditions  to  be  satisfied  for  such  applications  are  those  of  a 
clamped  edge.  The  clamped  edge  boundary  conditions  are  w{a)  =  0,  wtR{a)  =  0,  and  u(a)  =  0  or,  in  terms 
of  p,  25(1)  =  0,  u5,p(l)  =  0,  and  2(1)  =  0.  The  first  of  the  boundary  conditions  leads,  as  in  the  previous 
problem,  to  the  same  forms  of  the  solutions  given  in  (4.158)  and  (4.120).  The  coefficients  Cx  and  C2  in  this 
case,  however,  must  now  be  solutions  of  the  two  boundary  conditions  2(1)  =  0  and  w,p(l)  =  0,  which  can 
be  written  as 


CllCi  +  C12C2  —  (1  —  1^4  )^T 

Ciber'^x)  +  C^bei^fci)  =  0, 


(4.166) 

(4.167) 
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w(R)  -vs-  R  Comparison 

Pinned  (Hinged)  Edge  Boundary  Conditions  (  Sc  =  -5  KPa,  Ss  =  20  Pa) 


Figure  9:  Comparison  of  theory  to  geometrically  linear  FE  ^-displacement  results  (pinned  edge). 


where  (4.167)  follows  from  (4.104),  and  cn  =  J>n,  ci2  =  P12,  where  pn  and  p12  are  given  by  (4.162)  and 
(4.163),  respectively.  Equations  (4.166)  and  (4.167)  are  easily  solved,  yielding 

c,=  (1 

In  this  case,  both  the  radial  and  axial  displacement  solutions  contain  the  common  load  factor 

tl  = r  ■  (4-170) 

which  is  the  only  occurrence  of  the  intrinsic  stresses  in  the  solutions.  By  choosing  the  parameters  in  this 
factor  appropriately,  i.e.,  such  that 

tl  =  r  —  =  0,  or  equivalently,  No  =  J\f  —  —  +  ^  Q  =  0,  (4.171) 

2  K/q 

it  should  be  possible  to  achieve  a  state  of  no  deformation  from  the  initial  parabolic  shape.  For  example,  if 
there  is  no  pressure  difference  then  one  should  be  able  to  adjust  the  coating  stress  (or  one  of  the  coating 
stresses  of  a  multilayer  coating)  to  satisfy  A f  -  (70  9/ 2«o)  a  =  0.  In  a  Og  environment,  this  condition  reduces 
to  A f  =  0,  which  for  a  single  coating  is  equivalent  to  having  a  coating  stress  given  by  the  simple  prescription 

Sc  =  ~Ss.  (4.172) 

tlc 


bei'(Ari) 


Ciibei'(fci)  -  ci2ber,(fci)  J  ’ 

ber^fci) 


Lcnbei'(&i)  -ci2ber'(*i)J 


(4.168) 

(4.169) 
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Equation  (4.171),  and  its  special  case  (4.172),  define  what  we  refer  to  as  “on-design”  prescriptions  for  a 
coating  stress  that  will  maintain  the  initially  parabolic  shape  of  a  coated  membrane  after  removal  from  the 
mold  upon  which  it  was  cast  and  coated.  Note  that  this  simple  prescription  holds  only  for  clamped  edge 
boundary  conditions.  If  the  coating  stress  is  less  than  the  on-design  value,  we  say  that  the  membrane  is 
undercompensated.  On  the  other  hand,  if  the  coating  stress  is  greater  than  the  on-design  value,  we  say  that 
the  membrane  is  overcompensated. 

The  clamped-edge  solutions  are  perhaps  the  most  important  for  the  analysis  of  near  net-shape  coated 
membranes  used  as  optical  quality  reflectors.  For  this  reason,  we  have  carried  out  a  more  extensive  com¬ 
parison  with  finite  element  models  than  for  the  previous  two  boundary  value  problems.  The  details  of  this 
work  were  presented  at  the  43rd  AIAA  Structural,  Structural  Dynamics,  and  Materials  Conference  in  April 
2002  [20],  and  accepted  for  publication  in  2003  [21].  Here,  we  reproduce  comparisons  of  our  geometrically 
linear  theory  to  both  geometrically  linear  and  nonlinear  finite  element  results. 

The  model  in  this  case  is  a  10  m  diameter  (a  =  5  m),  //2  coated  membrane,  in  a  0 g  environment.  The 
material  and  geometrical  properties  of  the  coating  and  membrane  are  the  same  as  the  last  two  boundary 
value  problems,  viz.,  Ec  =  44  GPa,  Es  =  2.2  GPa,  uc  =  vs  =  0.4,  coating  thickness  hc  =  l  pm,  and 
membrane  thickness  hs  =  20  pm.  We  have  somewhat  arbitrarily  specified  a  membrane  CTE  mismatch  stress 
of  11  MPa,  which  corresponds  to  an  on-design  coating  stress  of  —220  MPa,  calculated  from  equation  (4.172). 

In  Figures  (10)-(13)  we  show  the  predicted  effects  on  the  displacement  components  of  overcompensating 
the  membrane  stress  by  either  1%,  shown  in  Figures  (10)  and  (11),  or  10%,  shown  in  Figures  (12)  and  (13). 
An  edge  effect  beginning  some  10  to  20  cm  from  the  edge  is  observed  in  both  the  theoretical  and  FE  results. 
The  graphs  of  radial  displacement  pass  through  the  origin  and  are  linear  until  the  onset  of  this  edge  effect, 
while  the  axial  displacement  curves  are  quite  flat  until  the  edge  effect  occurs.  To  enhance  the  detail  of  these 
edge  effects,  we  have  shown  in  the  Figures  only  the  final  meter  before  reaching  the  edge.  As  in  the  previous 
cases  considered,  the  agreement  between  theoretical  and  geometrically  linear  FE  results  is  excellent. 


Figure  10:  Comparison  of  theory  and  finite  element  results  for  radial  displacement  u(R)  when  the  coating 
stress  is  1%  off-design. 


However,  in  Figure  1 1  the  geometrically  nonlinear  axial  displacement  is  roughly  97%  of  the  theoretical 
prediction  for  the  1%  off-design  case,  while  in  Figure  13  the  nonlinear  prediction  is  only  72%  of  the  theoretical 
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w(R)  -vs-  R 

Intrinsic  Coating  Stress  1%  Off-Design 


Figure  11:  Comparison  of  theory  and  finite  element  results  for  axial  displacement  w(R)  when  the  coating 
stress  is  1%  off-design. 


one.  This  tendency  for  the  theoretical  result  (and  linear  FE  results)  to  overestimate  the  axial  displacement  is 
made  evident  in  Figure  14,  which  indicates  that  for  this  example  geometrical  nonlinearities  become  important 
when  the  coating  stress  is  more  than  2%  off-design,  while  the  theory  is  fairly  accurate  when  off-design  by 
less  than  2%. 


4.3  General  Solution  for  an  Initially  Flat  Laminate 

The  solution  procedure  of  §4.2  included  several  steps  involving  division  by  the  dimensionless  parameter  kq 
defined  by  (4.64).  For  an  initially  flat  laminate,  requiring  k  =  0  hence  k0  =  0,  these  divisions  cannot  be 
made.  The  solutions  for  a  flat  laminate  presumably  follow  from  those  of  the  parabolic  laminate  in  the  limit 
k0  -»•  0  (i.e.,  the  focal  length  /  ->  oo),  but  we  prefer  to  return  to  equations  (4.66)  and  (4.68),  prior  to  any 
divisions  by  kq,  and  solve  them  anew.  Setting  kq  =  0  in  these  two  equations  yields 


\xr  T  Xg  fix/  (/Sr  +  «*)]. 


and 

-A  (ttr  +  ),p  +  Q  2  =  °> 

each  of  which  can  be  immediately  integrated  to  obtain 

Xj~  4“  xg  fix/  (/Sf*  4”  /s$)  C]_ ,  (4.173) 

and 

q  p 

Kr  +  Kg  =  C2  4-  , 

where  C\  and  C2  are  arbitrary  integration  constants.  Using  (4.77),  equation  (4.174)  can  be  written  as 


(4.174) 


wD  q  p 

w,pp  H - -  C2  4-  x 


A  4 


(4.175) 
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Figure  12:  Comparison  of  theory  and  finite  element  results  for  radial  displacement  u(R)  when  the  coating; 
stress  is  10%  off-design.  ° 


Multiplying  (4.175)  through  by  p  yields 


Pw,pp  +  u5,p  =  (pu>,p)  —  C2p  +  “ 


.3 


*PP  ’  •“  ~  ^  r  «  ^4 

which  can  be  integrated  twice  to  obtain  the  general  solution  for  w(p): 

&(p)  =  c4  +  c3  lnp  +  c2 

4  A  64 


(4.176) 

where  c3  and  c4  are  arbitrary  integration  constants.  We  must  set  c3  =  0  in  order  for  the  solution  to  be 
regular  at  p  —  0,  hence  the  general  solution  regular  at  the  origin  is 

~(  \  P2  Q  P4 

™{p)=ci+c2-  +  (4.177) 

Returning  to  equation  (4.173),  we  use  equations  (5.39)  and  (5.40)  to  replace  xr  and  xe,  yielding 

xr  +  xg  =  y_Va  (€°rr  +  fee)  +  (b  +  K)  (nr  +  ne)  =  a  +  /?„  ( nr  +  ks)  , 

from  which 

crr  +  €ee  =  ci  (1  -  vA)  -  (1  -  v\)  b  (/cr  +  ne) ,  (4.178) 

where  to  used  the  definition  (4.50)  to  replace  ftv.  Equation  (4.174)  can  be  used  to  eliminate  (/cr  +  )  in 

(4.178),  yielding 

€°rr  +  4>e  =  ci  (1  -  vA)  -  (1  -  v\)  q  (4.179) 

Substituting  in  (4.179)  for  the  strain  components  from  (4.39)  (with  f\p  =  -n0p  =  0),  we  obtain  the 
following  differential  equation  for  u(p): 

2- p  +  “  =  ci  (J  "  va)  ~  (1  “  vX)  Q  £  J- 


38 


w(R)  -vs-  R 

Intrinsic  Coating  Stress  10%  Off-Design 


Figure  13:  Comparison  of  theory  and  finite  element  results  for  axial  displacement  w(R)  when  the  coating 
stress  is  10%  off-design. 


Multiplying  through  by  p  and  integrating  yields  the  general  solution 

u(p)  =  Cj  +  Cl  (1  -  VA)  P-  -  (1  -  ul)  q  A  ^  (4.180) 

where  C5  is  a  final  integration  constant.  We  must  set  C5  =  0  to  obtain  a  solution  regular  at  p  =  0: 

u(p)  =  Cl  (1  -  va)  I  -  (1  -  "a)  9  ^  fg-  (4.181) 

Restoring  the  dimensional  variables  and  constants  in  equations  (4.181)  and  (4.177),  we  write  the  general 
solutions  for  u(R)  and  w(R)  regular  at  the  origin  as 

u(R)  =  e°R-  quR3,  qu  =  (4‘182) 

and 

/  ir>\  k  r,2  .  r,4  _  A(P+1o9) 

w(R)  =■  w0  -  -  R  +  qwR  ,  qw  =  —  _  Biy  (4.183) 

where  we  have  introduced  new  arbitrary  constants  e°,  wo  and  fc,  as  well  as  coefficients  qu  and  qw  containing 
the  pressure  and  gravity  loads.  The  three  arbitrary  constants  e°,  wq  and  k  must  be  determined  by  applying 
boundary  conditions. 

4.3.1  Flat  Pressurized  Laminate,  Clamped  at  the  Edge 

For  a  flat  laminate  under  pressure,  the  boundary  conditions  insuring  that  the  edge  R  —  a  of  the  laminate 
remains  fixed  are  u(a)  =  0,  w(a)  =  0,  and  u/(a)  =  0.  Applying  the  first  of  these  to  (4.182)  yields  an 
expression  for  e°: 

e°  =  qua 2, 
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Coating  Stress  Percent  Off-Design  (%) 

Figure  14:  Comparison  of  theory  and  finite  element  results  for  apex  displacement  w(0)  as  a  function  of  the 
percent  that  coating  stress  is  off-design. 


hence  the  clamped-edge  solution  for  the  radial  displacement  is 


“<*>  -  -  r!)  =  -  fl2>' 


(4.184) 


The  first  boundary  condition  for  the  axial  displacement  yields  from  (4.183)  the  following  expression  for  the 
apex  displacement  wq: 


k 

wo  =  -a2  -  qwa4, 


which  brings  the  solution  to  the  form 


(4.185) 


w(R)  =  |  (a2  -  f?2)  -  qw  (a4  -  R4) . 
Setting  the  derivative  of  this  expression,  evaluated  at  R  =  a,  to  zero  yields 

w'(a)  =  -ka  +  4qwa3  =  0, 


from  which 


k  =  4  qw  a2 


Substituting  this  result  in  (4.185)  then  gives 


_ _  „4  _  A(p+~f0g)  4 

“  64  (AD-BO"' 

Using  k  from  (4.187)  in  (4.186),  the  clamped-edge  solution  for  the  axial  displacement  reduces  to 


(4.186) 


(4.187) 


(4.188) 


(4.189) 
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We  note  that  for  a  single  material,  the  coefficient  B  appearing  in  (4.189)  and  defined  in  equation  (1.3)  is 
zero.  The  axial  displacement  solution  in  that  case  reduces  to  the  one  given  in  equation  (62),  p.  55  of  [22]  (if 
the  effects  of  gravity  are  ignored).  The  general  effect  of  the  multilayers  is  to  reduce  the  composite  bending 
stiffness  D  by  an  amount  B2 /A.  It  is  also  worth  remarking  that  these  solutions  have  no  dependence  on  the 
residual  stress  loads  Si .  Thus,  in  the  absence  of  a  pressure  load,  and  if  gravity  is  ignored,  then  qu  -  0  and 
qw  =  0  and  both  displacements  vanish  for  clamped-edge  boundary  conditions  regardless  of  the  residual  stress 
levels. 

4.3.2  Unpressurized  Laminate  with  Free  Edge:  Generalized  Stoney  Formula 

We  suppose  now  that  there  is  no  pressure  load,  and  that  the  effects  of  gravity  can  be  neglected.  Then  qw  =  0 
and  qu  =  0  in  the  general  solutions  (4.182)  and  (4.183)  for  the  displacements,  reducing  them  to 

u(R)=e°R ,  (4.190) 

and 

w(R)  =  w0  -  ^R2.  (4.191) 

From  these  solutions  and  the  definitions  (4.18)  and  (4.19)  with  T,r  =  0,  we  find  immediately  that 

£rr  —  £©©  =  €°,  kRR  —  &©©  =  k,  (4.192) 

hence  the  stress  resultants  and  couples  (4.28)-(4.31)  reduce  to 

Nr  =  Nq  =  N  =  Af  ■+■  (-4  4-  Ay)  e°  4*  (B  4*  BJ)  k,  (4.193) 

and 

Mr  =  Me  =  M  =  -M  -  (B  +  Bv)  e°  -  (D  4-  Dv)  k,  (4.194) 

where  N  and  M  are  constants.  We  note  that  from  the  definitions  (1.2)-(1.6),  the  summed  coefficients 
appearing  in  (4.193)  and  (4.194)  take  the  forms 


N 

A  +  Av  =  ^  Mu 

i- 1 

(4.195) 

-  N- 1  N 

B  +  By  =  $  E  E  hihk  ^  ~  Bk )  > 

i  =  1  ~  i-f  1 

(4.196) 

1  N 

D  +  Dv  =  12  E  hiBi  [h2i  +  3  (h  -  -  12^-1  -  6)]  > 

i=l 

(4.197) 

where  we  have  introduced 

Tp ' 

Bi  =  Qi  (1  +  Vi)  = 

(4.198) 

and  made  use  of  the  first  definition  in  (1.8),  i.e., 

Q.  -  Bi 

Qt  1  -vV 

(4.199) 
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to  obtain  the  last  equality  of  (4.198).  As  mentioned  earlier  in  §3.3,  the  quantities  Qt  and  B,  are,  somewhat 
confusingly,  both  referred  to  in  the  literature  as  the  biaxial  modulus  of  the  material  in  layer  i.  Henceforth,  we 
reserve  the  name  biaxial  modulus  for  Bi,  defined  in  (4.198).  Note  that  the  ratio  of  biaxial  moduli  occurred 
earlier  in  the  condition  (4.155)  for  no  axial  deformation  of  a  free-edge  membrane  simply  supported  at  its 
center. 

We  now  consider  a  coated  membrane  laminate  with  free-edge  boundary  conditions  NR(a)  =  0  and 
Mr{o)  -  0.  These  two  conditions  allow  the  determination  of  the  arbitrary  constants  e°  and  k  appear¬ 
ing  in  the  general  solutions  (4.190)  and  (4.191).  In  fact,  we  have  found  in  equations  (4.193)  and  (4.194)  that 
both  NR  and  Mr  are  constants  for  this  problem,  so  that  c°  and  k  are  simply  the  algebraic  solutions  of  the 
linear  system  in  these  two  unknowns  obtained  by  setting  NR(a)  =  N  =  0  and  MR(a )  =  M  =  0: 

0  =  Af  +  {A  +  Av)  t°  +  (B  +  Bv )  k, 

0  =  -M  -  (B  +  Bv)  €°  -  (D  +  Du)  k. 

The  solutions  are  easily  found  to  be 

e°  =  -(.D  +  D*)M+{B  +  Bv)M 
{A  +  AU){D  +  DV)-{B  +  BU)2' 


(4.200) 

(4.201) 

(4.202) 


,  _  {B  +  By)  A f  —  (A  +  Ay)  Ad 

“  (A  +  AV)(D  +  D„)  -  (B  +  B„) 2'  (4-203) 

This  expression  for  k  represents  the  generalization  to  a  multilayer  coated  substrate,  with  arbitrary  layer 
thicknesses,  of  an  important  result  due  originally  to  G.  G.  Stoney  [23].  Note  that,  as  discussed  in  §1 
following  equation  (1.20),  k  is  the  negative  vertex  curvature  of  the  paraboloid  defined  in  this  case  by  (4.191). 
When  applied  to  a  substrate  with  a  single  coating,  it  reduces  to  a  result  discussed  in  [24,  pp.  140-143],  and 
can  also  be  found  (in  various  forms)  in  several  recent  publications  [25,  26,  27,  15,  28]. 

The  remaining  integration  constant  w0  =  w( 0)  is  the  axial  displacement  of  the  center  R  =  0.  We  assume 
that  there  is  at  least  one  value  of  R,  say  Rq,  at  which  the  coated  membrane  is  simply  supported,  meaning 
that  tn(.Ro)  =  0,  hence  w0  =  kRl/2.  For  support  at  the  center,  Rq  =  0,  we  have  w0  =  0,  and 


Wssc{R)  =  -2R2' 


simply  supported  at  the  center, 


(4.204) 


and  for  support  at  the  edge,  Rq  =  a,  we  have  w0  =  ka2 /2,  hence 


lOsse  (R)  —  —  (a2  —  R2) , 


simply  supported  at  the  edge. 


(4.205) 


In  either  case,  w  (R)  =  —kR  and  w"(R)  =  —k,  so  if  k  >  0  the  axial  displacement  will  have  its  maximum 
value  of  wsse( 0)  =  ka2 / 2  at  the  center,  while  for  k  <  0  it  will  have  its  minimum  value  of  it'sse(0)  =  0  at 
the  center  (and  the  positive  value  wssc(a)  =  —ka2 / 2  at  its  edge).  The  maximum  axial  deflection,  which  will 
occur  either  at  the  center  or  at  the  edge,  is  thus 


w  -  —  a2 

wmax  —  2  U  ■ 


(4.206) 


Note  that  under  these  boundary  conditions  the  requirements  for  no  displacement  of  the  initially  flat  coated 
membrane  are  e°  =  0  and  k  =  0,  which  imply  N  =  0  and  M  =  0,  the  same  conditions  found  in  (4.154)  for 
no  displacement  of  an  initially  parabolic  membrane  satisfying  the  same  boundary  conditions. 

These  results  can  be  applied  to  the  design  of  multilayer  coatings  that  will  compensate  for  curvature  in¬ 
duced  by  residual  stresses  in  the  coatings.  For  example,  Cao,  et  al  [29],  discuss  such  compensation  techniques 
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for  reducing  or  eliminating  unwanted  curvature  due  to  intrinsic  coating  stresses  in  high-reflectance  thin-film 
micromirrors,  and  Liu  and  Talghader  [30]  discuss  compensating  for  micromirror  curvature  due  to  thermal 
stresses  caused  by  CTE  mismatches  between  layers.  In  both  papers  an  extension  [31]  to  multilayer  stacks 
of  Timoshenko’s  [32]  one-dimensional  analysis  of  bimetal  thermostats  is  used  for  a  preliminary  analysis  of 
the  residual  stress  effects,  replacing  each  instance  of  a  modulus  E  by  the  biaxial  modulus  Ej (1  —  v)  (but 
continuing  to  denote  this  biaxial  modulus  by  the  same  symbol  E ).  For  comparison  with  Figure  3  of  Ref¬ 
erence  [29]  we  have  plotted  in  Figure  15,  below,  the  absolute  value  |fc|  of  our  vertex  curvature  (4.203)  as  a 
function  of  the  thickness  hi  (in  nanometers)  of  an  outer  dielectric  coating  of  Si02,  placed  over  a  gold  coating 
of  thickness  /i2  =  0.5  pm  on  a  polysilicon  membrane  substrate  of  thickness  h3  =  3.5  pm.  We  have  used  their 
data  for  the  various  layers,  viz.,  Ex  =  80GPa,  St  =  100  MPa  for  the  Si02  outer  coating,  E2  =  80GPa, 
52  =  -100  MPa  for  the  interior  gold  coating,  and  E3  =  125  GPa,  S3  =  0  for  the  polysilicon  substrate.  The 
curvature  corresponding  to  a  displacement  from  an  initially  flat  surface  of  wmax  — .  |A/10|  defines  a  “flatness 


M  -vs-V 


Figure  15:  Absolute  value  |fc|  as  a  function  of  thickness  of  outer  Si02  dielectric  coating  layer. 


window”  for  the  allowable  thickness  variation  of  the  Si02  coating.  This  curvature  can  be  found  from  (4.206), 
i.e., 

W  =  2^  =  5^.  (-*-207) 

which  requires  knowledge  of  the  specimen  size.  In  Reference  [29]  the  specimens  were  square  with  sides  of 
length  200  pm,  corresponding  roughly  to  circular  specimens  of  radius  a  =  100  pm.  The  curvature  corre¬ 
sponding  to  this  radius  is 


(4.208) 
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so  that  for  visible  wavelengths  on  the  order  of  500  nm,  we  have 

1*1  «  10,  (4.209) 

corresponding  to  the  line  |fc|  =  10  in  Figure  15.  We  observe  that  for  this  particular  case  our  flatness  window 
is  about  63  nm,  compared  to  the  much  narrower  20  nm  flatness  window^  showm  in  Figure  3  of  [29].  Based  on 
correspondence  with  Professor  Talghader,  it  is  believed  that  the  discrepancy  is  due  to  either  a  typographical 
error  in  the  data  reported  for  Figure  3,  or  an  error  in  the  simulator  being  used  at  the  time  of  publication.  A 
rerun  of  their  current  simulator  with  the  data  used  here  shows  excellent  agreement  with  the  result  reported 
in  our  Figure  15.  It  should  be  noted  that  the  thickness  required  to  flatten  the  mirror  is  generally  not  the 
same  as  that  which  would  zero  out  the  radial  expansion  or  contraction  given  by  (4.190)  and  (4.202).  In 
Figure  16  w'e  show  the  variation  of  radial  edge  deflection  u(a)  =  e°a  in  nanometers  as  the  Si02  thickness 
increases,  using  the  same  data  as  that  of  Figure  15.  The  flatness  window  corresponds  to  a  maximum  radial 
expansion  of  about  3  nm,  and  a  minimum  of  1.5  nm. 


Figure  16:  Radial  edge  displacement  u(a)  as  a  function  of  thickness  of  outer  Si02  dielectric  coating  layer. 


5  Geometrically  Nonlinear  Shell  Laminate  Problems 

In  this  final  Section  we  present  the  geometrically  nonlinear  theory  of  a  coated  membrane  shell.  Such  prob¬ 
lems  are  notoriously  resistant  to  analysis  in  terms  of  known  elementary  or  special  functions.  We  specialize 
immediately  to  the  axially  symmetric  and  initially  flat  subset  of  these  problems.  Our  hope  was  to  generalize 
to  stress-coated  multilayered  plates  some  of  the  earlier  work  [33,  34]  on  plates  of  a  single  material  with  no 
residual  stress.  In  the  work  cited,  Chia  [33,  pp.  112-115]  applied  perturbation  theory  to  find  approximate 
solutions  of  the  plate  problem,  w'hile  Way  [34]  attacked  the  problem  with  power  series  expansions,  a  tech¬ 
nique  that  had  been  used  much  earlier  by  Hencky  [7]  in  solving  the  equations  of  membrane  theory.  As  we 
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shall  see,  the  existence  of  residual  stresses  and  the  appearance  of  additional  terms  due  to  the  multilayer 
structure,  greatly  complicate  the  solution  by  either  perturbation  series  or  power  series  methods.  Our  efforts 
to  this  point  have  been  for  the  most  part  unsuccessful,  and  we  must  be  content  to  simply  record  what  has 
been  accomplished  in  the  hope  that  it  may  be  useful  as  a  reference  point  for  others  interested  in  these  types 
of  problems. 

As  in  §4,  the  leading  order  displacement  components  for  a  nonlinear  shell  laminate  have  the  Kirchhoff- 
Love  forms  given  in  equation  (8.31.1)  of  Reference  [1]: 

Uz  —  w,  Ur  =  u  -  Zw,R,  U@  =  v  -  Z  (5.1) 

where  it,  v,  and  w  are  functions  of  R  and  0  only.  The  strain  components  given  in  equations  (8.35.1)— (8.40.1), 
however,  contain  nonlinear  terms  in  the  displacement  component  derivatives: 
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where  the  Z-independent  terms  of  the  last  three  equations  are  given  by 
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The  equilibrium  equations  are  given  in  terms  of  stress  resultants  by  equations  (8.130.I)-(8.132.I): 

( RNr)r  —  Nq  +  NRete  —  0, 


(5.4) 

(5.5) 

(5.6) 

(5.7) 

(5.8) 


(R2NRe)R  +  RN@,e  =  0,  (5.9) 

and 

| R  (oj,rNr  +  +  Oit)]  R  +  (u>,rNrq  +  — ©  +  Qe)  q  +  (p  +  7o  <?)  R  =  0,  (5.10) 

where  we  recall  that  u(R)  =  w(R,Q)  +  T(R).  The  shear  stress  resultants  Qr  and  Q©  appearing  in  the 
last  equilibrium  equation  can  be  eliminated  in  favor  of  stress  couples,  which  yields  the  alternative  equation 
(8.140.1)  of  Reference  [1]: 

\r  (u,rNr  +  ^NRe)  +  {RMr)r  -  Me  +  Mfle.e]  R 

r  i  j  (5.ii) 

+  Iw^IVr©  +  — j^-iV©  +  (R2  Mr©)  r  +  —  +  (p  +  ■yo  g)  R  =  0. 
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The  stress  resultants  and  couples  appearing  in  the  equilibrium  equations  are  given  by  equations  18  141 II- 
(8.146.1):  v  •  •  ; 


Nr  =  N  +  AeRR  +  4„ee0  +  BkRR  +  B„kee,  (5.12) 

Nrq  —  AeeRe  +  J5©  fcfl© ,  (5.13) 

Nq  =  N  +  A„e°RR  +  A  e0e  +  B„kRR  +  Bk&e,  (5.14) 

Mr  =  -M  -  BeRR  -  Bv€%q  -  DkRR  -  Dvkee,  (5.15) 

Mrq  =  -Be  e°RQ  -  DekRe,  (5.16) 

Me  =  —M  —  Bve°RR  —  5e00  -  DukRR  -  Dkee ,  (5-17) 


w  hich  are  the  same  as  those  of  the  previous  Section,  except  that  the  strains  contain  nonlinear  terms,  according 
to  equations  (5.5)-(5.7).  The  constants  appearing  here,  as  elsewhere,  are  defined  in  §1. 

5.1  Reduction  to  an  Axisymmetric  System 

The  reduction  to  an  axisymmetric  system  follows  along  the  same  lines  as  §4.1,  i.e.,  we  assume  that  all 
variables  are  independent  of  ©,  and  find  from  equilibrium  equation  (5.9)  that 

Nrq  =  0,  for  all  R,  (5.18) 

hence  from  (5.13),  since  kRe  =  0  from  (5.5)  (and  .4©  is,  in  general,  nonzero): 

ffi©  =  2  [V’R  ~  r)  =  (519) 

The  last  equation  can  be  integrated  to  give  v{R)  =  v0R,  where  vq  is  an  arbitrary  constant,  hence  if  v  vanishes 
for  any  non-zero  value  of  R,  then  v  vanishes  for  all  R.  These  results  imply  from  (5.16)  that 

Mrb  =  0,  for  all  R.  (5.20) 

Noting  that  equation  (5.11)  can  be  integrated  once  (and  the  integration  constant  discarded  to  insure  regu¬ 
larity),  the  system  of  equations  to  be  solved  thereby  reduces  to  the  following: 

{RNr)  r  -  Ne  =  0, 

( RMr)  r  —  Me  +  R  (wtR  +  r*)  Nr  +  (p  +  703)  —■  =  0. 

where 

Nr  =  J\f  +  AeRR  +  A^tQQ  +  B  kRR  +  B^kee, 

Ne  =  Af  +  A„eRR  +  Ae@Q  +  BrkRR  +  Bkee, 

Mr  =  -M  -  BeRR  -  Bv e00  -  DkRR  -  D„kee, 


(5.21) 

(5.22) 

(5.23) 

(5.24) 

(5.25) 
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.M©  —  — A4  —  Bv  €^jrj  —  B  600  —  Dvkjm  Z^fc©©, 


and 


<4*  =  u,R  +  r  j*w  *  +  \tfR ,  foil  =  w,m, 


(5.26) 

(5.27) 


0  _  U  r  _  ^ 

C©0  “  ^  * 


(5.28) 


As  in  §4.1,  we  transform  these  equations  to  dimensionless  forms  by  introducing  a  dimensionless  coor¬ 
dinate  p,  dimensionless  displacement  components  u  and  w ,  and  a  dimensionless  initially  curved  reference 
surface  T,  defined  by 

p  =  “ ,  U  =  -,  w  =  — ,  r  =  -,  (5.29) 

~  a  a 


together  with  the  following  dimensionless  constants: 


^  = 


A- 


JV 

T  “  Eft’ 


M  =  w, 


A’  -  -  h  v-  '  -  1ft’  ■  r!  - 

and  four  dimensionless  dependent  variables  xr,  xg,  yT,  and  j/#  defined  by 


9  = 


(p  +Jog)a 
Eh 


Xr  — 


Nr-  AT 


_  Ne-Af 


_  Mr  +  M 

Wha  ’ 


Eh  ’  "  Eh 

We  also  introduce  dimensionless  curvatures  /tr  and  Kg  defined  by 


M©  +  Af 
Ve  =  ~^T — ■ 
Eha 


Kr  =  altRR  =  w, 


w . 


,pp 


Kg  =  ak@@  —  — — , 
P 


and  four  new  dimensionless  constants 

B 


b  = 


Eha 


K  s 


By 

Eha’ 


d  =  =- 


£> 


£/ia2 


£ha2’ 


(5.30) 


(5.31) 


(5.32) 


(5.33) 


noting  that  the  strain  components  are  already  dimensionless,  and  have  the  following  forms  in  terms  of  the 
dimensionless  variables  defined  above: 


efii?  =  £©e  -  “• 

The  strain  components  must  satisfy  a  modified  form  of  the  compatibility  condition  (3.40),  viz., 

Peee,p  ~  €rr  “  ee©  ~  F,p  wtP  -  -w*p, 

and  there  is  an  additional  compatibility  condition  involving  the  dimensionless  curvatures: 

P^6yp  =  ^ r  R'Bi 


(5.34) 


(5.35) 


(5.36) 


which  is  easily  derived  from  the  definitions  (5.32). 

Substituting  these  dimensionless  quantities  into  the  equilibrium  equations  (5.21)  and  (5.22),  we  reduce 
them  to 

xr,p  +  -  (xr  -  xg )  =  0,  (5.37) 


and 


l lr,p  +  -  (  Vr  -  ye  )  +  (w,p  +  f  >J5)  (r  +  Xr)  +.  9 1  =  o, 


(5.38) 
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respectively.  The  constitutive  relations  (5.23)-(5.26)  take  the  forms 


Xr 


(e°RR  +  Coe)  +  bnr  +  bvKg, 


(5.39) 


1 

l-ul 


(va£°rr  +  e©©)  +  bv  Kr  4-  6/cp, 


(5.40) 


Vr  =  -be°RR  -  b„e%e  -  Pkt  -  d„Kg,  (5.41) 

Do  =  ~b^eRR  —  5e©©  —  dv  kt  —  dug.  (5-42) 

Equations  (5.39)  and  (5.40)  are  easily  solved  for  e<RR  and  e©e  in  terms  of  the  x’s  and  k’s: 

t°RR  =  Xr  -  vAxg  -  Pkt  -  Pvk6,  (5.43) 

e©0  =  xe  -  vAxr  -  Pvkt  -fine,  (5.44) 

where  we  recall  the  definitions  first  given  in  (4.50): 

P  =  b  -  vAb„,  pv  =  £>„  -  uAb.  (5.45) 

We  note  here  that  the  dimensionless  radial  displacement  component  can  be  expressed  in  terms  of  xr  and  xe , 
and  Kr  and  Kg  (which  themselves  depend  only  on  derivatives  of  w) ,  by  combining  (5.34)  and  (5.44)  to  obtain 

u(p)  =  p(xg  -  vA  xT  -  pv  Kr  -  P  k0)  .  (5.46) 


Taking  the  derivative  with  respect  to  p  of  (5.44)  and  substituting  the  result,  together  with  (5.43)  and 
(5.44),  in  the  c-compatibility  relation  (5.35),  we  obtain 

P  {x0,p  —  vAxr,p  —  Pv  Krtp  —  Pne,p)  =  (1  +  VA){xr  —  Xg)  —  (/?  —  Pv)  {kt  —  Kg)  —  TtP  Wtp  —  ^V?p- 
In  this  expression,  we  use  (5.36)  to  replace  (kt  —  /eg),  and  (5.37)  to  replace  (xr  —  xg),  which  yields  the 


following  result  for  the  compatibility  condition: 

p  [xr  +  xg  -  pv  ( Kr  +  Kg)]>fi  =  -f  ,p  -  -V?p.  (5.47) 

Next,  we  substitute  (5.43)  and  (5.44)  in  (5.41)  and  (5.42)  to  obtain 

yr  =  -pxr  -  puxg  -  6  Kr  -  6vKg,  (5.48) 

yg  =  -PvXr  -  Pxg  -  SuKr  -  6  Kg,  (5.49) 

where  we  introduced  another  pair  of  constants: 

&  —  d  —  bp  —  bypt, ,  8V  =  du  bp„  —  bvp.  (5.50) 

The  derivative  of  yr  with  respect  to  p  yields  from  (5.48): 

Vr,p  =  ~PxT,p  Pvx6,p  ~  $  Kr,p  ~  bvKg^p.  (5.51) 
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Erom  the  last  three  equations  we  obtain 


Vr,p  p  ( 2/r 


i)  —  ft  %r,p  fiv%6,p  &  ftr,p 


-  i[(/?  -  Pv){xT  -  Xg)  -  (6  -  S„)  (Kr  -  He)]- 


In  this  expression,  we  again  use  (5.36)  to  replace  («r  -  kq ),  and  (5.37)  to  replace  (xr  —  x$),  which  simplifies 
it  to 

Ur,p  4“  “  ( Ur  Ho'}  =  fiv  {%r  “I"  ^  S  (kt  -j-  Kq  )^  .  (5.52) 

Substitution  of  (5.52)  into  (5.38)  yields  the  following  form  of  the  axial  equilibrium  equation: 

-Pv  (xr  +  xe)iP  -  6  (Kr  +  Ke)tP  +  (w,p  +  f,p)  (r  +  xr)  +  q  |  =  0.  (5.53) 


Equation  (5.37)  can  be  used  to  eliminate  x#  in  (5.47)  and  (5.53),  and  the  problem  we  are  left  with  is  to 
solve  the  coupled  differential  equations  (5.47)  and  (5.53)  for  w  and  xr.  The  radial  displacement  u  is  then 
determined  from  (5.46). 


5.2  Approximate  Solutions  for  Initially  Flat  Laminates  Using  Perturbation 
Methods 

Here,  we  apply  perturbation  methods  to  an  initially  flat  coated  membrane  laminate,  for  which  T(p)  =  0 .  The 
problems  we  consider  will  be  of  two  distinct  types.  In  the  first  we  are  concerned  with  finding  approximate 
pressure  -vs-  axial  displacement  relations  for  fixed  residual  stress  loads  defined  by  r.  In  the  second,  the 
pressure  and  gravity  loads  defined  by  q  are  assumed  absent,  and  we  investigate  the  compressive  intrinsic 
coating  stress  contained  in  r  required  to  produce  buckling,  i.e.,  out-of-plane  deflections  of  the  laminate  from 
its  initially  flat  state. 


5.2.1  Pressure  Versus  Axial  Displacement  Curves 

For  the  first  type  of  problem,  we  follow  Chia  [33],  assuming  series  expansions  for  the  dimensionless  variables 
(which,  however,  are  not  the  same  as  Chia’s  dimensionless  variables)  in  a  perturbation  parameter  oj  =  w( 0), 
the  central  axial  displacement,  viz., 

w(p)  =  wn(p)cjn  =  ™i(p)Ct;  +  ™2(p)u2  +  ws(p)uj3  + -  (5.54) 

n= 1 


Xr  (p)  =  Y  xn(p)  0Jn  =  Xi(p)  U>  +  X2  (p)  U2  +  X3(p)  Ul3  +  . . . , 
n—l 

and 

q  =  ^2  QnUn  =  qi  u  +  q3  w2  +  q3  w3  +  . . . . 

n—  1 

Denoting  a  derivative  with  respect  to  p  by  a  “prime” ,  we  then  have 

=  Y  =  Y  Xr’p =  Y 

n~  1  n=l  n~  1 

The  nonlinear  terms  appearing  in  equations  (5.47)  and  (5.53)  are  given  to  fourth  order  in  u  by 
(w,p)2  =  (w[)2w2  +  2w[  w'2  w3  +  [{w'2)2  +  2u>i  ?4]  u>4  +  . . . , 


(5.55) 


(5.56) 


(5.57) 


(5.58) 
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and 


(w,p)  far  +  t)  —  w[  r  cj  +  [w[  (x\  +  r)  +  w*2  r)]  u2  +  [iyj  (£2  +  r)  +  tyj  (xi  +  r)  +  «4  r]  u;3 

+  K  (*3  +  t)  +  w2  (x2  +r) +  w'3(x!  +t) +  w'4r]  u4  +  ...,  (5.59) 

respectively.  When  these  perturbation  series  are  substituted  in  (5.47)  and  (5.53),  we  obtain  the  following 
equations  as  coefficients  of  the  indicated  powers  of  oj  for  the  first  four  orders: 


w‘:  < 


u2:  < 


'  p(pA  +  2*0'  -  pvp  (w'{  +  f)'  =  0, 

.  Pv  ( Pxi  +  )#  +  S  (w"  +  —  i  gjp  —  w[  r  =  0, 

C  p  (i Px2  +  2i2)'  -  &p  (<  +  f)'  +  1  (wj)2  =  o, 

[  Pv  (px 2  +  2x2)  +  S  (w2  +  —  \  q2p  -  w[  (xj  +  t)  -  w'2t  =  0, 


(5.60) 


(5.61) 


L) 


3  . 


P  (p4  +  2x3y  -  p„p  (ic"  +  f)'  +  w[w'2  =  0, 

< 

.  $v  (px  3  +  2  x3)f  +  S  (w3  +  -  \q3p  -  w[  (x2  +  t)  -  w'2  (xi  +  r)  -  w'3t  =  0, 


(5.62) 


P  ( Px4  +  2*4)'  -  Pi,  p  (w'l  +  |  [(^2)2  +  2u>(  iCjj  =  0, 

u>4  :  4  , 

Pv  (px 4  +  2x4)  +  2  |  94  p  —  u;[  (X3  +  r) 

,  —  w2  (x2  +  r)  -  w'3  (xi  +  t)  —  w'4t  =  0. 

The  boundary  conditions  for  a  rigidly  clamped  edge  at  R  =  a,  or  p  =  1,  are 

2(1)  =  0,  23(1)  =  0,  and  23iP(l)  =  0. 

Recalling  equation  (5.46)  for  u(p)  in  terms  of  xr,  xg,  and  23,  and  using  xg  =  pxTtP  +  xr  from  (5.37)  to 
replace  xg  in  (5.46),  the  first  boundary  condition  of  (5.64)  is  equivalent  to 

XrA !)  +  (1  “  *a)  ®r (1)  -  P„W,pp{l)  =  0,  (5.65) 

where  we  used  the  third  boundary  condition  to  eliminate  the  term  proportional  to  p.  Substitution  of  the 
perturbation  series  (5.54)-(5.57)  in  these  boundary  conditions  yields  the  rigidly  clamped  edge  boundary 
conditions  on  the  perturbation  coefficient  functions,  viz., 

<(1)  +  (1  -  vA)  xn(l)  -  &,<(!)  =  0,  u»n(l)  =  0,  and  <(1)  =  0.  (5.66) 

In  addition,  we  have  the  so-called  “false”  [35,  p.  167]  boundary  conditions  on  the  functions  wn  evaluated  on 
axis  at  p  =  0: 


(5.63) 


(5.64) 


u>i(0)  =  1,  and  in„(0)  =  0,  for  1,  (5.67) 

which  follow  from  the  requirement  that  the  apex  displacement  25(0)  =  ui  in  (5.54). 
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As  our  first  application  (although  somewhat  of  a  digression),  we  consider  the  specialization  of  these 
perturbation  equations  to  true  membranes.  Comparing  equations  (3.43)  and  (3.39)  of  §3.2  to  equations 
(5.47)  and  (5.53),  respectively  (with  f(p)  =  0),  we  see  that  the  geometrically  nonlinear  membrane  theory 
follows  from  the  latter  equations  by  setting  <5  =  0  and  (3„  =  0.  The  system  (5.60)-(5.63)  of  perturbation 
equations  thereby  reduce  to 


p  (px[  +  2x1)'  =  0, 
-  |  q1p  -  w[  t  =  0, 


(5.68) 


p  ( px'2  +  2<e2)'  +  |  (wj)2  =  0, 
-\q2P  -  (xi  +  t)  -  tt>2 r  =  0, 


(5.69) 


p  (px'3  +  2x3)'  +  w[w2  =  0, 

-  |  qzp  -  w'i  (x2  +  t)  -  w'2  (xi  +  t)  -  w'3  t  =  0, 


(5.70) 


p  (px 4  +  2x4 )  +  2 


(w'2)2  +  2w'1w'3 


=  0, 


[  -  I  q4p  -  w[  (x3  +  r)  -  w'2  (x2  +  r)  -  w3  ( xi  +  r)  -  w'4  r  =  0, 


(5.71) 


p  ( px'5  +  2x5)'  +  (w[  w4  +  w'2w3)  =  0, 

-\q3p  -  w[  (x4  +  r)  -  w'2  (x3  +  r)  -  w'3  (x2  +  r)  -  w\  (xr  +  r)  -  w'5  r  =  0, 


(5.72) 


where  we  have  appended  the  fifth-order  equations.  The  boundary  conditions  (5.66)  for  a  rigidly  clamped 
edge  reduce  to 


i^(l)  +  (1  -  pa)  xn(l)  =  0,  w„(l)  =  0, 


(5.73) 


noting  that  the  boundary  condition  ui,p(l)  =  0  of  (5.64)  does  not  apply  to  a  true  membrane.  The  apex 
displacement  conditions  (5.67)  are  unchanged.  In  the  process  of  solving  these  equations  for  the  functions 
Xi(p),  the  following  identity  is  useful: 


(px'i  +  2  Xi)' 


(5.74) 


Beginning  with  the  first  of  equations  (5.68),  two  integrations  give  the  general  solution 

Xi  (p)  =  cx  +  ^,  (5.75) 

where  c\  and  C2  are  arbitrary  integration  constants.  We  must  set  C2  =  0  to  insure  that  x\  is  regular  at  p  =  0. 
The  first  boundary  condition  of  (5.73),  with  n  =  1,  then  yields  ci  =  0  hence 

Xi  (p)  =  0,  Vp.  (5.76) 


The  second  of  equations  (5.68)  can  be  integrated  to  give  the  general  solution 

wi(p)  =  ~f^p2  +  °3, 


(5.77) 
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(5.78) 


and  the  first  condition  of  (5.67)  yields  C3  =  1,  hence 


wi(p)  =  1  ~  f~P2’ 

4  r 

The  second  boundary  condition  of  (5.73),  with  n  =  1,  then  yields 

0i  =  4t,  (5.79) 

and  the  w\  -solution  (5.78)  reduces  to 

(p)  =  1  -  p2,  =>  w[  =  -2 p.  (5.80) 

Substituting  for  w[  in  the  first  of  equations  (5.69),  two  integrations  give  the  general  solution 

*a(/>)  =  ci  “  \p2i  =►  *2  =  -£/>■  (5.81) 

where  one  integration  constant  has  been  set  to  zero  to  insure  regularity.  The  first  boundary  condition  of 
(5.73),  with  n  =  2,  then  yields 

1  (S-vA\ 

C,  =  4(j^)’  (5  82) 

which  can  be  substituted  in  (5.81)  to  obtain 

=  J  (l^  -  "O '  <5  83) 

Substituting  for  X\  and  w[  in  the  second  of  equations  (5.69),  and  integrating  the  result  once,  yields 

W2(p)  =  (l  -  |^)p2  +  c2,  (5.84) 

and  the  condition  (5.67)  with  n  =  2  forces  c2  =  0.  The  final  boundary  condition  u;2  ( 1 )  =  0  then  gives 

02  =  4r,  (5.85) 

which  yields  from  (5.84): 

w2(p)  =  0,  Vp.  (5.86) 

Continuing,  we  substitute  w'2  =  0  from  (5.86)  in  the  first  of  equations  (5.69),  which  leads  to  the  solution 

x3(p)  =  0,  Vp,  (5.87) 

after  applying  the  first  boundary  condition  of  (5.73).  Substituting  for  xu  x2,  w[  and  w'2  in  the  second  of 
equations  (5.69)  yields  an  equation  which  integrates  to 


(5.88) 


after  setting  the  integration  constant  to  zero  to  satisfy  u>3(0)  =  0.  The  boundary  condition  uj3(1)  =  0  then 
gives 


93  =  4r  + 


1  /5  -vA\ 

2  VI  -vA)' 


(5.89) 
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which  yields  from  (5.88): 


=  ^y (>-/>’).  '=*  =  i-v)- 


Proceeding  in  the  same  way  for  the  next  two  orders  we  obtain  the  following  solutions: 


i  r  i  +  v  a 


48  r  [l  —  i/ a 


-  P2  ( 3  —  2  p2)  , 


94  =  4t, 


for  n  =  4,  and 


t(p)  =  (!  -  P2)  >  =*  w4  =  (1  ~  2P2)’ 


q5  =  4r  + 


1  / 19  +  5  i/a 


144  7  \  \  —  V A 


1  f 61— 251/a \  ,2  ,  1  /^9-5t/a\  4 

^5(P)  576 r2  l  1  —  i/a  J  P  +  64 r2  (  1  -  vA  )  P 


144  t2 


for  n  =  5.  Note  that  w^{p)  =  -w3{p),  and  x5{p)  =  -x4(p).  Substituting  the  solutions  qit  i  =  1, . . . ,5  into 
the  perturbation  expansion  (5.56)  for  the  dimensionless  load  q,  we  obtain 


1  (5 -va 


q  =  4rw(l+u  +  u2+w3  +  w4)  +  -  (  j 


"s  +  m-  «s. 

144  t  \  1  —  va  ) 


The  sum  occurring  in  the  first  term  on  the  right-hand  side  is  a  finite  geometric  series,  hence  we  can  write 
(5.97)  in  the  form 


( l-tu5\  1  (5 -va\  3  1  / 19  +  5  i/a  \  5 

5  =  4r"  (t^t)  +  2  (rr^j  "  +  144?  "  ■ 


(5.98) 


Restoring  dimensions  using  the  definitions  in  equations  (5.29)  and  (5.30),  noting  that  to  =  wo/a ,  we 
write  this  in  terms  of  the  physical  variables  (ignoring  the  effects  of  gravity)  as 


h  Wo  \J\f  (1  —  (Jl 


IV  OL/y  I  J. 

p  =  4 - —  I  — 

a  a  h  \  1 


1  (5  -  i/a^  p  (w0\2  1  E2  h  ( 19  +  5  va^  (Wo}4 

\  a  )  ~  576 


8  \1  — i/a/  Va/  576  N  \  1  —  i/a 


(?)‘  ■ 


Equation  (5.99)  is  similar  to  Equation  (3)  of  the  Beams  paper  [8]  cited  in  §3.3,  if  we  replace  his  To  by  J\T/h) 
and  replace  our  finite  geometric  series  by  1.  However,  the  coefficients  of  his  third  and  fifth-order  terms  are  a 
factor  of  2  greater  than  ours,  and  there  is  a  difference  in  the  dependence  on  Poisson’s  ratio  in  the  fifth-order 
term. 

Returning  now  to  the  general  plate  perturbation  equations  (5.60)-(5.63),  we  assume  that  the  tensile 
residual  membrane  stress  dominates  r,  so  that  r  is  positive.  It  follows  that  the  parameter  72  defined  by 


where  A  =  (5  +  /?2  =  d  -  (l  -  v\)  62, 


(5.100) 
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is  also  positive  (the  constant  A  was  introduced  earlier  in  equation  (4.69)).  It  is  convenient,  then,  to  introduce 
a  new  independent  variable 


*  =  7  P>  (5.101) 

and  new'  perturbation  functions  defined  by 

Wnfa)  =  Wn(x/y),  Xn(x)  =  Xn(x/7).  (5.102) 

Our  system  of  perturbation  equations  can  be  written  in  terms  of  these  new  functions,  and  the  variable  x,  as 

|  x  (xX[  +  2  X1)'  -  pvy2x  ( W?  +  *£)'  =  0, 

(xX[  +  2X 0'  +  <572  (w?  +  ^)'  -  ftx  -  W[  t  =  0, 


or  : 


(5.103) 


a;2: 


+  2  X2)'  -  P„y2x  (w?  +  +  172  {W[f  =  0, 

/?„  (xX2  +  2X2)  +  J72  (w£  +  ^  x  —  W[  (Xi  +  r)  —  W2  t  =  0, 


f  x  (xlj  +  2X3)'  -  Pvy2x  (w"  +  K)'  +  72  WJ  =  0, 


(5.104) 


(xX3  +  2X3)'  +  <572  (w"  +  *?)'  - 

-  W'  (X2  +  r)  -  W2'  (Xi  +  r)  -  =  0, 


I  x  (iXJ  +  2X4)'  -  fay2  x  (Wf  +  ^)'  +  |72  [(^  )2  +  2  W[W^  =  0, 


(5.105) 


w4:  < 


(5.106) 


&  (xX'  +  2X4)'  +  6 y2  (W”  +  J&)'  -  ^  x  -  W{  (X3  +  r) 

-  WJ,  (x2  +  t)  -  (Xj  +  r)  -  W[  r  =  0. 

The  boundary  conditions  for  a  rigidly  clamped  edge  at  p  =  1,  i.e,  x  =  -y,  are 

7X;(7)  +  (1  -  pA)  X„(7)  -  72&  W"(7)  =  0,  Wn{y)  =  0,  and  W'n{y)  =  0,  (5.107) 

and  we  must  also  satisfy  the  false  boundary  conditions  on  the  functions  Wn  evaluated  on  axis  at  p  =  0, 

equivalent  to  x  —  0: 


Wi{ 0)  =  1,  and  Wn{ 0)  =  0,  for  n  #  1. 
From  the  first  equation  of  (5.103)  we  have 

(xX(  +  2X0'  =  Pv 72  (Wiw  +  ^)\ 
This  can  be  substituted  in  the  second  equation  of  (5.103)  to  obtain 


(5.108) 


(5.109) 
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(5.110) 


where  we  used  (5.100).  This  equation  can  be  integrated  once,  and  written  as 

Wl 

W{'  +  — i  -  Wl  =  Qi  X2  +  C3, 

X 


where  C3  is  an  arbitrary  integration  constant,  and 


gi 

7  2r  ’ 


(5.111) 


The  solution  regular  at  x  =  0  of  the  homogeneous  differential  equation  associated  with  (5.110)  is  Iq(x ), 
where  Iq  is  the  modified  Bessel  function  of  order  zero,  i.e., 


Wlh(x)  -  Cxloix),  (5.112) 

where  Ci  is  another  integration  constant.  A  particular  solution  is  easily  found  to  be 

Wlp(x)  =  -QlX2  +  C2,  (5.113) 

where  C2  is  a  new  arbitrary  constant.  The  general  solution  is  the  sum  of  the  homogeneous  and  particular 
solutions: 


Wi(x)  =  Ciloix)  -  Qix2  +  C2.  (5.114) 

The  conditions  to  be  satisfied  by  W\  are  obtained  by  setting  n  =  1  in  (5.107)  and  (5.108).  Thus, 

Wi(7)  =  0,  W{(  7)  =  0,  and  Wi(0)  =  1.  (5.115) 

The  third  condition  yields  C2  =  1  —  C\ ,  and  the  first  condition  then  gives 

1  -  <3i72 


Cx  = 


Since  C2  =  1  —  C\ ,  we  thus  have 


C2  = 


1  -iohY 

Qa2-io(i) 


(5.116) 


(5.117) 


i-ioh) 

The  final  boundary  condition,  W{ (7)  =  0,  will  determine  Qi  or,  equivalently,  q\  .  From  (5.114)  we  have 

W[{x)  =  Cih(x)  -  2QlX,  (5.118) 

noting  the  identity 


Evaluating  at  x  =  7  yields 

0  =  CMi)  -  2Qi7  = 

which  can  be  solved  for  Qi.  Thus, 

Qi  = 

and 


Io(x)  =  h(x). 

1  -  <3172 


1  -  io(7) 


h( 7)  “  2<2i7, 


■^1(7) 


7  [7 ^1(7)  +  2  -  2  7o(7)]  72  L  9  72  1296 


10  16  7 2  4\ 

+  77  +  7^7  +  °(T  ) 


qi  =  4  72t 


hh) 


ilbhil)  +  2  -  2/0(7)]] 


=  645  +  ^-t  +  t  0(72), 

y 


(5.119) 

(5.120) 

(5.121) 

(5.122) 
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where  the  second  equality  made  use  of  the  following  power  series  representations  for  the  Bessel  functions:. 


J0(7)  =  i  +  I72  +  J_74  +  _J_  7e  , 

1  4  7  64  7  64  •  36 7  ’ 

7 1\  (7)  =  -  72  +  —  74  -| - - —  76  + 

'  ;  2  7  16  7  64-6  7 


(5.123) 

(5.124) 


We  next  integrate  (5.109)  once  to  obtain 


xx[  +  2x,  =  /?„72  (w;'  +  ypj  +  C4, 


(5.125) 


where  C4  is  an  arbitrary  integration  constant.  From  (5.118)  we  obtain  for  the  second  derivative  of  Wx: 


using  the  identitity 


W"(z)  =  C,  /„(*)  -  t-hix)  -  2 Q1 


I[{x)  =  I0(x)  -  -h(x). 

x 


(5.126) 


(5.127) 


Substituting  from  (5.118)  and  (5.126)  into  (5.125),  we  find 

xX[  +  2X4  =  o(x)  +  C5, 

where  C5  is  a  new  constant.  Multiplying  thru  by  x  and  making  use  of  the  identitity 

xlo(x)  =  [z/^x)]', 

we  obtain  after  one  integration  the  general  solution  for  Xi(x): 

Xx{x)  =  /3,72C,^  +  %  +  % 
x  2  x2 

We  set  C6  =  0  to  insure  regularity  at  a;  =  0,  obtaining 

Xi{x)  =  &  7 2C1^  +  % 

x  2 


(5.128) 


(5.129) 


The  arbitrary  constant  C5  must  be  determined  by  the  first  boundary  condition  of  (5.107)  with  n  =  1,  viz., 


7^(7)  +  (1  -  va)  *1(7)  -  72^W/'(7)  =  0. 


From  (5.129)  we  have 


x[(x)  =  72  C\ 


Io(x)  _  2h(x) 
X  x2 


(5.130) 


(5.131) 


and  we  substitute  from  (5.131),  (5.129),  and  (5.126),  evaluated  at  x  =  7,  into  (5.130)  to  obtain  the  following 
expression  for  C5 : 


C5  =  2 


(5.132) 


Substitution  of  this  expression  into  (5.129)  completes  the  solution  of  the  first-order  perturbation  equations. 
In  the  second-order  equations,  we  have  from  the  first  equation  of  (5.104): 


{xX'2  +  2X2)'  =  &72  (V"  +  HV 


f)  ~ 


(5.133) 
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which  we  substitute  into  the  second  equation  of  (5.104),  obtaining 


(w"  +  ^'  -W^  =  2Q2x  +  Wl  +  (5.134) 


where 


Q,  =  i-f. 
4  7  2t 


(5.135) 


Substituting  in  the  right-hand  side  of  (5.134)  for  W[  and  Xx  from  (5.118)  and  (5.129),  respectively,  yields 
(w?  +  ^  -Wl  =  2  (Qr  -  Q2  -  +  hfQVj  x 


+  Ci 


Using  (5.119)  and  the  identity 


[J02(*)  —/?(*)]  s  2 


1  +  Yr  (WrfQi  ~  C6) 


>  _  JfW 


/i(s)  + 


3/?„72Ci2  J?(*) 


2r 


x 


(5.136) 


the  last  equation  can  be  integrated  once  to  obtain 


W' 

Wo'  +  —  -  w2  =  C6  + 


(Qi  - 


Q2  - 


CsQi  ,  /?, 


2  T 


+ 


V^Ql)  x2 


+  Ci  1  + 


2r 


(2^72Qi  -  C5 ) 


/o(£)  + 


3/?!/72Cl  r  _2 


4r 


[/2(x)  -  /2(x)] ,  (5.137) 


where  C6  is  an  arbitrary  integration  constant.  The  solution  of  the  homogeneous  equation  associated  with 

(5.137)  is  again  I0(x).  However,  we  have  been  unsuccessful  in  finding  an  appropriate  particular  solution,  due 
to  the  appearance  of  the  quadratic  products  of  Bessel  functions  in  the  last  term  on  the  right-hand  side  of 

(5.137) .  This  brings  to  a  rather  abrupt  end  our  attempt  to  solve  the  general  perturbation  equations  past  the 
first  order.  The  one  useful  result  is  perhaps  the  first-order  expression  (5.122)  for  the  dimensionless  pressure 
load: 


qi  =  4  72t 


h  (7) 


=  64(5  +  4-  rO(72), 


(5.138) 


1 7  [7  ^1(7)  +  2  -  2/0(7)]  J 

where  the  last  equality  was  obtained  by  using  power  series  representations  for  the  modified  Bessel  functions. 
Prom  the  definitions  of  5  and  r  we  have 


S  =  = 


Eha2 


D 


(AB2  -  2 AVBBV  +  ABl 


V 


A2  -  A2 


)]• 


T  Eh' 


(5.139) 


Substituting  (5.138)  in  (5.56)  yields  the  first-order  expression  for  pressure  as  a  function  of  the  dimensionless 
apex  displacement,  i.e.,  q  =  q\w,  which  can  be  written  using  (5.139)  and  the  definitions  of  q  and  w  (to 
restore  the  dimensions)  as 


pa 

Th 


f  64  r  _  (AB2  -  2AVBBV  +  ABl  ^ 
{Eha2  [V 


V 


A2  -  A2 


J  J 


40  M  2n 

+  JTh+T°h) 


1  ^ 
J  a  ' 


(5.140) 


When  specialized  to  a  single  material  in  which  there  is  no  residual  stress  (hence  B  =  0,  Bv  =  0,  and  AT  =  0 ), 
equation  (5.140)  reduces  to  a  form  of  Chia’s  first-order  result: 


pa  _  (  64  D  \  wo 
Eh,  ~  \Eha2 )  V’ 


(5.141) 
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noting  that  D  —  Eh3 /[12{\  - 1/2)]  for  a  single  material.  This  also  follows  directly  from  the  linear  plate  theory 
of  a  single  material,  as  evidenced  by  evaluating  equation  (4.189)  on-axis  at  R  =  0. 

It  should  be  noted  that  one  of  the  major  sources  of  our  difficulty  involves  the  appearance  of  terms 
containing  the  coefficient  /?„.  We  have  shown  in  equation  (4.60),  and  its  specialization  (4.61)  to  two  layers, 
that  this  coefficient  can  be  written  as  a  sum  of  terms,  each  of  which  depends  upon  a  difference  of  two  Poisson’s 
ratios.  Thus,  if  all  layers  have  the  same  Poisson  ratio,  this  coefficient  must  vanish,  considerably  simplifying 
the  resulting  equations.  It  is  unlikely  that  this  will  ever  be  exactly  true,  but  it  is  perhaps  worthwhile 
investigating  the  perturbation  solutions  obtained  by  making  the  assumption  that  terms  containing  /3„  can  be 
set  equal  to  zero.  We  have  found  that  with  such  an  assumption  we  gain  a  solution  through  second  order,  but 
are  unable  to  solve  the  third-order  equations,  and  so  it  is  here  that  we  end  our  discussion  of  the  perturbation 
series  approach  to  finding  approximate  pressure  -vs-  axial  displacement  curves. 

5.2.2  Buckling  Due  to  Compressive  Intrinsic  Stress  Loads 

For  the  second  type  of  perturbation  problem,  we  assume  the  pressure  and  gravity  loads  defined  by  q  to  be 
absent,  and  instead  investigate  the  compressive  intrinsic  coating  stress  contained  in  r  required  to  produce 
out-of-plane  deflections  of  the  laminate  from  its  initially  flat  state.  The  perturbation  series  expansions  for 
w(p)  and  xr{p)  are  the  same  as  before,  i.e.,  from  (5.54)  and  (5.55): 

=  YlWn^u}n  =  wi(p)  w  +  w2(p)  w2  -I-  w3{p)  ui3  +  . . . .  (5.142) 

n=l 


xr(p)  -  X)z„(p)wn  =  x1{p)oj  +  x2(p)u2 +x3(p)u>3  +  ...,  (5.143) 

n=l 

but  we  now  set  q  =  0  in  the  equilibrium  equations,  and  construct  a  perturbation  expansion  for  the  intrinsic 
stress  load  r ,  assumed  to  be  negative,  beginning  with  the  n  =  0  term: 

T  =  Tnun  =  ~T0  -  ri  w  -  r2  ui2  -  r3  u>3  - -  (5.144) 

n=0 

The  clamped  edge  boundary  conditions  on  the  perturbation  functions,  as  well  as  the  false  boundary  condi¬ 
tions,  are  also  the  same  as  in  (5.66)  and  (5.67),  i.e., 

xn(!)  +  (!  ~  vA)  a;n(l)  -  &w"(l)  =  0,  wn(  1)  =  0,  and  wj,(l)  =  0,  (5.145) 

wi(0)  =  1,  and  tu„(0)  =  0,  for  n  ^  1.  (5.146) 

The  nonlinear  term  appearing  in  equation  (5.47)  is  given  to  fourth  order  in  w,  as  earlier,  by 

(w,p)2  =  {w\ )2  w2  +  2  w[  w2  w3  +  [(u4)2  +  2w[  Wg]  u>4  +  —  (5.147) 

but  the  nonlinear  term  in  (5.53)  now  involves  the  perturbation  series  for  r,  and  takes  the  form 

{wtP)  ( xr  +  t)  =  -w[  T0  w  +  [u;J  (ii  -  n)  -  u>2  r0]  w2+ 

[W1  (x2  ~  T2)  +  w'2  (Xi  -  Ti)  -  w'3  T0]  U3 

+  K  (2:3  -  r3)  +  w'2  ( x2  -  t2)  +  w'3  (X!  -  Ti)  -  w\  To]  u4  +  .  .  .  .  (5.148) 


When  these  perturbation  series  are  substituted  in  (5.47)  and  (5.53),  we  obtain  the  following  equations  as 
coefficients  of  the  indicated  powers  of  u  for  the  first  four  orders: 


ui 


1  . 


’  P  {px[  +  2x0'  -  Pvp  (<  +  f  )'  =  o, 

i 

k  Pv  {px i  +  2x0'  +  6  (w"  +  +  w[  To  =  0, 


(5.149) 
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P  (px'2  +  2  x2)'  -  Pvp  (w2  +  +  \  ( wi )2  =  °> 

OJ2  :  < 

0V  (px'2  +  2x2)'  +  s  (w2  +  -  w[  (xi  -  Ti)  +  w'2  To  =  0, 

p  (px3  +  2 x3)'  -  fivP  [w3  +  ^)  +  w'iW'2  =  0, 

0V  ( px'3  +  2x3)'  +  6  (w%  +  ^)  -  w[  ( x2  -  T2)  -  w'2  (xi  -  Ti)  +  w'3t0  =  0, 

p  ( px\  +  2x4)'  -  Pvp  (w4  +  +  \  ( w'2f  +  2 w[w'3  =  0, 

Pv  (px 4  +  2x4)'  +  6  [w'l  +  ^)  -  w[  (x3  -  r3) 

k  -W2  (x2  -  T2)  -  w3  (xi  -  n)  +  w'4t0  =  0. 


(5.150) 


(5.151) 


(5.152) 


Beginning  with  equations  (5.149),  we  note  that  the  first  can  be  used  to  eliminate  X\  in  the  second,  yielding 
the  following  equation  for  w\  after  integrating  once: 


//  w  T0  ^ 
U)"  +  —  +  “T*  wi  =  Co, 
P  & 


(5.153) 


where  Co  is  an  arbitrary  integration  constant,  and  A  was  introduced  earlier  in  equation  (4.69),  repeated 
here: 

A  =  J  +  0t  =  d  -  (1  -  v\)  b\  (5.154) 

Setting 

A2  =  g,  (5.155) 


equation  (5.153)  can  be  written  as 


w"  +  —  +  A2wi  =  Co- 
P 


(5.156) 


The  solution  regular  at  the  origin  of  the  homogeneous  equation  associated  with  (5.156)  is  Jo(Ap),  the  ordinary 
Bessel  function  of  order  0,  and  the  particular  solution  is  simply  Co/A2  =  C2,  hence  the  general  solution  is 


wi{p)  =  Ci  Jo(Ap)  +  C2  • 


(5.157) 


This  must  satisfy  the  boundary  conditions  (5.145)  and  (5.146)  with  n  =  1.  The  conditions  Wi(0)  =  1  and 
1*4(1)  =  0  yield  the  solutions 


c‘  =  r rrj^ry  = 

The  derivative  of  (5.157)  gives 

w[  (p)  =  -CiXJi(Xp), 

and  the  boundary  condition  tnj(l)  =  0  then  requires  A  to  be  a  solution  of 

Ji(A)  =  0, 


(5.158) 


(5.159) 


(5.160) 


that  is,  A  must  be  a  zero  of  Jj,  the  Bessel  function  of  first  order.  The  first  (non-zero)  zero  of  J\  occurs  at 
A  =  3.832,  hence  the  smallest  value  of  tq  is,  from  (5.155): 


T0  =  A2  A  =  14.684  A. 


(5.161) 
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Thus,  from  (5.161)  and  (5.144),  the  critical  value  of  r  above  which  there  will  be  an  out-of-plane  displacement 
is  ! 


r  =  -r0  =  -14.684  A.  (5.162) 

Restoring  the  dimensional  constants  using  the  definitions  in  §5.1,  equation  (5.162)  yields  the  following  value 
for  the  critical  intrinsic  stress  resultant: 


AT0  =  - 


14.684 

a2 


(5.163) 


For  a  single-layer  coating,  we  have  A/"o  —  hcSc o  +  haSs,  hence  the  critical  coating  stress  Sc o  required  for  an 
out-of-plane  displacement  is 


14.684 
hca 2 


(5.164) 


Although  we  have  succeeded  in  continuing  the  solutions  of  the  perturbation  series  to  higher  orders,  it  is 
this  result  for  the  critical  coating  stress  that  is  of  greatest  interest  to  us.  We  thus  end  our  discussion  of 
perturbation  methods  here,  and  proceed  to  a  brief  discussion  of  solutions  in  the  form  of  power  series. 


5.3  Power  Series  Solutions 

In  1934,  Way  [34]  published  “exact”  power  series  solutions  for  the  geometrically  nonlinear  plate.  In  order  to 
follow  his  work  as  closely  as  possible,  we  introduce  his  normalizations  of  the  radial  variable  and  displacement 
components  with  respect  to  the  thickness  h  rather  than  the  radius  o,  denoting  them  by 


f?  a  tt  -  u  a  ~  u/  w  a  - 

{=  h~k U=h  =  Hu’  Wsk=j,w- 


(5.165) 

iSote  the  correspondence  with  Way’s  notation:  our  f  is  his  u-coordinate  and  our  radial  displacement  u  is  his 
variable  p.  Our  variables  xT  and  x$  correspond  directly  (in  the  absence  of  intrinsic  stress  loads)  to  Way’s 
variables  Sr  and  St ,  respectively,  which  we  distinguish  here  by  Xr  and  respectively,  so  that 


Xr  =  * =  sir- 


V  _  _  Ws-V 

=  ~WJT- 


We  also  introduce  new  dimensionless  curvatures  Kr  and  Kg  defined  by 


Kr  =  hkR  R  =  W'X  =  -Kr  Kg  =  hkee  =  ^  =  -ks. 


(5.166) 


(5.167) 


Way  did  not  consider  the  effects  of  gravity.  His  dimensionless  pressure  load,  which  we  denote  by  Q  (not  to 
be  confused  with  Way’s  scale-invariant  pressure  load  Q)  is  related  to  our  q  (in  the  absence  of  gravitational 
effects)  by 


o  =  L  =  - 

W  ~  E  aq’ 


(5.168) 

The  dimensionless  stress  load  r  does  not  appear  in  his  equations,  since  he  is  not  concerned  with  intrinsic 
stresses. 

The  equations  we  must  solve  are  the  compatibility  equation  (5.47)  and  axial  equilibrium  equation  (5.53) 
viz.,  5 


p  [xr  +  Xg  -  (Kr  +  Kg)]  p  = 

( xr  +  x$)  p  —  6  (/cr  +  p  +  w,p(t  +  av)  +  q  r  =  0, 


and 


(5.169) 

(5.170) 
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using  (5.37)  to  eliminate  xe ,  i.e., 

xe  =  (pxr)tP.  (5.171) 

In  terms  of  Way’s  variables,  these  take  the  forms 

*  [Xr  +  Xe-Bv  (Kr  +  Ke)]4  =  -\w\.  (5.172) 

and 

-Bv  (Xr  +  Xt)4  -V(Kr  +  Ke  )4  +  WA  (r  +  Xr)  +  Q  f  =  0,  (5.173) 

where 

X,  =  (£Xr)}C,=  Xr  +£Xr,c  (5.174) 

and  we  have  introduced  new  constants 


using  definitions  given  earlier.  Note  that  the  new  constants  do  not  depend  on  the  radius  a,  and  that  the 
terms  proportional  to  Bv  do  not  appear  in  Way’s  equations,  as  Bv  and  B  are  zero  for  a  plate  consisting  of  a 
single  material.  Also,  for  a  single  material  D  reduces  to  the  usual  bending  stiffness  Eh?/[  12(1  -  i/2)],  hence 
V  reduces  to  1/[12(1  -  v2)].  It  is  easy  to  show  that  the  radial  displacement  (5.46)  in  Way’s  variables  takes 
the  form 

U  =  £  {Xe  ~  uA  Xr  -  Bv  Kr  -  BKe ) ,  (5.176) 

where  the  constant  B  is  defined  by 

B  =  (5.177) 

Following  Way  [34],  we  seek  power  series  solutions  for  Xr  and  having  the  forms 

oo 

Xr(0  =  -r+^b2nen,  (5.178) 

n=0 

oo 

WM  =  Vs  e2n+1-  (5.179) 

n— 0 

The  series  for  W  is  found  by  integrating  (5.179)  term-by-term  to  yield 

W(0  =  Co  +  Vs  (— ^)  C2n+1  ^2n+2’  (5-180) 

where  cq  is  an  arbitrary  integration  constant  (the  value  of  W  when  £  =  0).  To  determine  the  coefficients  &2n 
and  C2n+i?  n  =  0, 1, 2 . . .,  we  note  that  from  (5.174)  and  (5.178)  we  have 

(oo  \  00 

-£r  +  b2nVn+1  =-t+  X)(2«  +  l )b2nen,  (5.181) 

n- 0  /  £  n= 0 

hence 

OO  1  OO 

(Xr  +  Xe  )4  =  -2r  +  2  ^  (n  +  1)  b2n  Vn  =  4  £  <n  +  x)  &2«  (5.182) 

n= 0  J  n— 0 

From  (5.179)  we  obtain 

oo 

WM  =  Vs  52  (2n  +  !)  C2n+1  Vn,  (5.183) 

71=0 
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hence  from  the  definitions  (5.167): 


n= o 


Kr  ”  -  ^8  E  (2n  +  1)  c2n+i  f2n,  ^0  =  =  V8  E  C2n+1  £2n. 


n=0 


From  these  two  expressions  it  follows  that 


(ffr  +  Kb)*  = 


E  (n  +  1)  c2n+1  f2n 


L  n=0 


=  4%/8  n(n  + l)c2n+i^2n  1. 


(5.184) 


(5.185) 


n=0 


Substituting  the  appropriate  series  expressions  in  the  compatibility  condition  (5.172)  and  axial  equilibrium 
equation  (5.173),  we  obtain 

00  .  /  PC  \  2 

E  <n  +  !)  ( ' b2n  ~  Vs  B„  c2n+1  U2n  +  £  c2n+1  C2n+1  =  0,  (5.186) 

n=0  V n=0  / 


and 


.  00  00 

-8  E  +  1)  (  Bv  b2n  +  V8  Pc2n+1  J  C2n  +  2  V8  E  C2"+i  ^"+1  *  E  f2"  +  0  f  =  0, 

"=°  n=0  - n 


(5.187) 


n=0 


respectively.  Using  Mathematica  again,  we  determined  the  following  recursion  relations  for  all  coefficients  in 
terms  of  the  two  unknown  coefficients  b0  and  Ci : 


.  _  v  Vs  Bv  „ 

2k+2  £>2k- 1-2  “1“  .  ~rrr\  02fc+3?  &  =  0,1,2,... 


C2fc+3  =  - 


(P  +  B2) 

VSBv 

8  (P  +  B2) 


P  +  B2) 


B2k+i  + 


(®  +  B2) 


C2*+3.  =  0, 1,2,  . . . 


(5.188) 

(5.189) 


where 


1  /  Q  \  i  *-i 

°3  =  8  [ui  +  boCl )  ’  C2fc+1  =  Pfc  +  1)2  _  !]  E  B2mC2k-y-2m,  k  =  2,3,4,...  (5.190) 


and 


2/c+2  =  — 


(*  +  !)(* +  2) 


E  ^2^1+1  ^2*1+1-2771,  fc  =  0,1,2,... 


m— 0 


(5.191) 


are  essentially  Way’s  [34]  coefficients,  with  his  factor  12(1  —  i^2)  omitted  in  the  numerators  of  the  C-coefficients 
(as  mentioned  earlier,  this  factor  corresponds  roughly  to  1/D  in  our  coefficients).  Our  unknown  coefficients 
correspond  directly  to  his  Bq  and  Ci,  i.e.,  bo  =  Bq  and  cx  =  Cx.  Thus,  for  example,  when  k  =  0  in  (5.188) 
and  (5.189)  we  obtain 


,  _  P  1  2  ,  y/8  Bv  1  f  Q  \ 

2  (P  +  S2)  2Cl  +  (P  +  B2)  8  (2V8  +  6°ClJ’ 

VW)  8  (wf  +  60  Cl)’ 


VS  B„  1  2 

^3  .  .in  \  .  Cl  *f 


8 (P  +  B3)  2  1  T  (P  +  B2)  8  V2V8 


(5.192) 

(5.193) 


where  we  used  (5.191)  with  k  =  0,  i.e.,  B2  =  -cf/2,  noting  that  c\  =  Cf. 
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The  coefficients  b0  and  C\  must  be  determined,  in  general,  from  the  boundary  conditions.  We  consider 
here  a  clamped  boundary,  for  which  the  boundary  conditions  at  the  circular  edge  R  =  a,  corresponding  to 
£  =  a/h  =  £0  (Way’s  constant  u0),  are 

!7(£o)  =  0,  W(£ o)  =  0,  and  Wc(£o)  =  0.  (5.194) 

Prom  equation  (5.176),  we  have 

u(t o)  =  6  [**(&)  -  VAXr(t o)  -  ft,  *.-(&)],  (5.195) 

noting  that  ifgi  is  proportional  to  W which  must  vanish  at  £  =  £o  according  to  the  second  boundary 
condition  of  (5.194).  Thus,  the  first  boundary  condition  corresponds  to 

Xe{£o)  -  ***r(&)  -  B,^r(£o)  =  0.  (5.196) 

Substituting  in  (5.196)  the  appropriate  series  expansions  from  (5.181),  (5.178),  and  (5.184),  evaluated  at 
£  =  £0,  yields  for  the  first  boundary  condition: 

o° 

£  [(2n  +  l-uA)b2n  -  n/8  Bv  (2n  +  1)  c2n+i  J  -  (l-i 'a)  r  =  0.  (5.197) 

n=0 

From  (5.179)  evaluated  at  £  =  £0,  the  third  boundary  condition  of  (5.194)  takes  the  simple  form 

f;  c2n+1e0n+i  =  o.  (5.198) 

71=0 

The  two  nonlinear  equations  (5.197)  and  (5.198)  must  be  solved  simultaneously  for  the  unknown  coefficients 
b0  and  c1;  and  then  the  second  boundary  condition  of  (5.194)  yields  the  remaining  unknown  coefficient  Co, 

i-e-, 

co  =  -V8±  (2^2)  C2n+ie0n+2,  (5.199) 

according  to  (5.180).  We  note  that  the  coefficients  bo  and  C\  are,  according  to  (5.178)  and  (5.183),  related 
to  the  values  of  Xr  and  on-axis  at  £  =  0  by 


bo  =  xr{  0)  +  r,  and  a  =  —  W^(0),  (5.200) 

respectively.  Since  is  proportional  to  the  curvature  of  the  axial  displacement,  we  expect  that  if  the 
pressure  load  changes  sign  (reversing  the  curvature),  so  also  must  Ci. 

The  solution  of  the  system  of  two  nonlinear  equations  (5.197)  and  (5.198)  for  b0  and  c\  is  a  daunting 
task,  unless  one  is  able  to  provide  good  initial  estimates  for  them.  Way  partially  avoided  this  difficulty  by 
first  assigning  values  to  bo  and  ci,  then  solving  (5.198)  for  the  corresponding  value  of  £o  =  a/h  (again,  his 
it0).  In  this  way,  the  second  boundary  condition  was  automatically  satisfied.  This  value  of  £o  must,  however, 
produce  a  radius  a,  for  a  given  specimen  thickness  h ,  near  the  actual  radii  of  test  articles  in  any  experiment 
to  which  the  theory  is  applied.  The  burden  of  solving  (5.197)  and  (5.198)  directly  is  thus  shifted  to  that  of 
choosing  values  for  b0  and  ci  that  lead  to  a  value  of  £o  which  provides  a  suitable  a  for  a  given  h.  The  solution 
of  this  difficulty  was  unfortunately  not  elaborated  upon  in  Reference  [34].  We  have  found,  however,  that 
for  values  of  £o  of  interest  to  Way,  i.e.,  for  £o  <  100,  approximate  values  for  bo  and  c\  can  be  obtained  by 
keeping  only  the  first  two  terms  of  the  boundary  conditions  (5.197)  and  (5.198),  and  solving  the  truncated 
equations  for  b0  and  c\  (this  procedure  leads  to  a  cubic  equation  for  Ci).  At  any  rate,  for  a  given  bo,  Way 
proceeded  to  determine  two  or  three  more  values  of  Ci  that  would  produce  values  of  £o  yielding  radii  in 
an  interval  containing  the  desired  radius  a.  Finally,  he  used  these  values  to  determine  by  interpolation  the 
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value  of  d  that  would  give  a  &  satisfying  the  first  boundary  condition  (5.197),  and  fixing  the  radius  a.  All 
other  quantities  could  be  determined  in  the  form  of  power  series,  using  the  values  for  b0  and  cj  satisfying 
the  required  boundary  conditions. 

Unfortunately,  values  of  £o  that  correspond  to  typical  coated  membrane  geometries  are  much  higher  than 
any  considered  by  Way.  For  example,  the  membrane  model  discussed  in  §3.3  has  a  radius  of  a  =  0.0762  m 
and  total  thickness  h  =  21  pm.  corresponding  to  £0  «  3629,  nearly  two  orders  of  magnitude  higher  than  the 
values  of  interest  to  Way.  For  large  aperture  coated  membrane  reflectors  having  radii  on  the  order  of  a  meter 
or  more,  fo  is  another  one  to  two  orders  of  magnitude  greater  still!  The  method  of  truncating  each  of  the 
two  boundary  conditions  at  the  second  term  and  solving  for  b0  and  cx  yields  completely  erroneous  estimates 
of  these  coefficients  for  such  large  values  of  £o-  We  have  been  thus  far  unsuccessful  in  developing  a  method 
for  estimating  and  ci  that  is  applicable  to  coated  membrane  geometries. 

5.3.1  Scale-Invariant  Functions  and  Constants 

Consider  a  transformation  of  the  dimensionless  coordinate  f  by  an  arbitrary  scale  factor  s,  generating  a 
scaled  coordinate  £: 

7  =  st-  (5.201) 

We  assume  that  the  functions  and  constants  appearing  in  equations  (5.172)-(5.174)  transform  to  new  func¬ 
tions  under  this  scaling  according  to  the  following  general  rules: 

W(0  =  s”W(0,  Kr(0  =  S1'  Kr(0,  KeCO  =  s'2  Ke(0,  (5.202) 

Xr(l)  =  s'*'  Xr(£),  Xg(£)  =  Sk*  Xg(S),  (5.203) 

B„  =  sm  Bu,  V  =  snX>,  T  =  sT  t,  Q  =  stQ,  (5.204) 

where  the  scale  factor  exponents  are  arbitrary.  For  any  scaled  function  /(£)  =  sa  /(£),  we  find  using 

the  chain  rule  for  differentiation  the  following  relations  between  derivatives  with  respect  to  the  scaled  and 
unsealed  coordinates: 

7,7  =  *  /,{  = 

Under  these  scalings,  then,  the  fundamental  equations  (5.172)-(5.174)  transform  to 

s~ll  [(^^r  +  sl-k*X0)  -  s-mB„  Kr  +  s^Ke)]^  =  -ls2~2pw\, 


-s~mBv  {sl~k'XT  +  sl~k'Xe)ri  -  s~n V  {sl~e'Kr  +  s1"'’]?*)- 

+  (s-rr  +  s~k'Xr)  +  s-'-'qI  =  0, 

2 

and 

sk>Xe  =  sk'Xr  +  sk'SXrl 

These  equations  can  be  simplified  to 

?  [{Xr  +  Sk'~k*Xg)  -  Bv  ( Sk'~rn-t'KT  +  Sk'-m-'*Kg)]-  = 


Bv  (Xr  +  Sk'  k*Xg)j  -  V  (s*i+"»-n-<j  Kr  +  Skl+m-n-(*  Kg)  - 

+  Wj  (skl+m~p-rT  +  Sm~pXr)  +  s*i+m-2-<g£  _  0) 


(5.205) 

(5.206) 

(5.207) 

(5.208) 

(5.209) 

(5.210) 
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and 


X,  =  *fcl-fc2(£Xr)if,  (5.211) 

respectively.  These  equations  will  retain  the  same  form  in  the  scaled  variables  as  in  the  unsealed  variables 
if,  and  only  if,  the  exponent  of  each  occurrence  of  the  scale  factor  s  is  zero.  The  equations  are  then  said  to 

be  form-invariant  under  the  scaling  transformation.  This  requirement  leads  to  a  linear  system  of  equations 

which  can  be  solved  for  eight  of  the  nine  exponents  in  terms  of  the  remaining  one.  For  example,  we  find  the 
following  solutions  in  terms  of  m :  fa  =  k2  =  r  =  2m  -  2;  lx  -  £2  =  m  -  2;  n  =  2m;  p  —  m;  and  t  =  3m  -  4. 
Thus,  the  fundamental  equations  are  form-invariant  under  the  following  one-parameter  group  (parameter  m) 
of  scaling  transformations: 

w®  =  *mW(0,  Kr(0  =  sm-2KT(0,  Keif)  =  sm~2Ke(0,  (5-212) 

Xr{f  =  s2m~2Xr(0,  Xeif  =  s2m-2Xe(0,  (5.213) 

t  =  s2m~2  t,  Bv  =  smBv,  V  =  $2mV,  Q  =  s3m~4Q.  (5.214) 

The  subgroup  of  this  transformation  that  preserves  the  material  constants  Bv  and  £>,  i.e.,  leaves  them 
scale-invariant,  is  obtained  by  setting  m  —  0.  In  this  case  the  dimensionless  axial  displacement  W(f)  is 
scale-invariant  as  well,  i.e.,  __  _  _ 

W(  0  =  W(  0,  Bv  =  Bv,  V  =  V.  (5.215) 

The  remaining  functions  and  constants  transform  under  this  subgroup  to 


Kr(0  =  s-2  Kr(S),  KeiS)  =  S~2  KeiO,  (5.216) 

XriO  =  s-2Xr(0,  XeiO  =  s~2Xe(0,  (5-217) 

t  —  s~2  t,  Q  =  s"4  Q.  (5.218) 

Using  (5.201)  to  eliminate  the  scale  factor  s  =  f/£  in  this  system,  we  find  the  following  scale-invariant 

combinations:  _ _  _ 

S2Kr(Z)  =  eKrit),  ?  Kt(t)  =  e  Kelt),  (5.219) 

t  Xrif)  =  exriO,  fXeCO  =  exe(0,  (5.220) 

fr  =  er ,  fQ  =  eQ.  (5.221) 


These  expressions  apply  at  all  points  including  the  edge  £o  =  cl/ h  hence,  in  particular,  Q£q  is  a  scale-invariant 
constant  (under  the  transformations  that  leave  Bv  and  V  scale-invariant) ,  and  the  stress-related  functions 
£2  XT  and  f2  X&  are  also  scale-invariant.  Way  [34]  used  Qt£  (his  qu g)  as  the  abscissa  in  several  graphs,  and 
discussed  scaling  properties  briefly  in  the  Appendix  of  [34]. 


6  Conclusions 

This  work  is  the  sequel  to  a  Report  [1]  in  which  asymptotic  methods  were  used  to  derive  theories  that 
would  aid  in  understanding  the  mechanical  behavior  of  a  stress-coated  membrane.  We  have  applied  those 
theories  to  a  number  of  boundary  value  problems,  obtaining  generalizations  of  well-known  solutions  for 
a  membrane,  plate,  or  shell  of  a  single  material  to  solutions  for  the  same  structure,  but  now  consisting 
of  a  multilayer-coated  polymer  material.  In  particular,  we  (i)  found  geometrically  linear  non-axisymmetric 
solutions  for  a  pressurized  membrane  laminate  having  a  general  ©-dependent  boundary  curve,  (ii)  determined 
linear  membrane  vibration  solutions  for  a  membrane  laminate  clamped  along  a  planar,  circular  boundary, 
(iii)  generalized  the  Hencky-Campbell  [6,  7]  geometrically  nonlinear  power  series  solutions  for  a  pressurized 
membrane  to  a  membrane  laminate,  and  discussed  the  application  of  these  solutions  to  the  interpretation  of 
bulge  test  data,  (iv)  derived  geometrically  linear  general  solutions  for  both  initially  parabolic  and  initially 
fiat  coated  membrane  laminates,  generalizing  the  Stoney  formula  for  such  a  structure  simply  supported  at 
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its  center,  and  (v)  initiated  a  search  for  perturbation  series  and  power  series  solutions  of  the  geometrically 
nonlinear  theory  of  an  initially  flat  coated  membrane  laminate. 

Perhaps  the  most  significant  accomplishment  of  this  work  was  the  discovery  of  simple  prescriptions  for 
the  coating  stress  that  would  maintain  the  shape  of  an  initially  parabolic  coated  membrane  after  removal 
from  the  mold  upon  which  it  was  cast.  The  details  of  the  solutions,  from  which  these  prescriptions  follow, 
are  presented  in  §4.2  and  in  particular  in  §4.2.4.  These  solutions  involve  linear  combinations  of  Kelvin 
functions.  Coating  prescriptions  are  given  there  for  membrane  laminates  both  with,  and  without,  pressure 
and  gravitational  loads.  The  prescriptions  are  presently  being  used  in  the  preliminary  design  of  a  near  net- 
shape  stress-coated  membrane,  and  it  is  hoped  that  such  a  structure  will  be  demonstrated  for  the  first  time 
in  the  very  near  future. 
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A  Elementary  Analysis  of  Stress  and  Strain  Due  to  CTE  Mis¬ 
match  Between  Membrane  and  Mandrel 


We  consider  a  polymeric  membrane  cast  on  a  mandrel,  and  allowed  to  cure.  The  membrane/mandrel  laminate 
is  then  subjected  to  an  annealing  process  in  which  it  is  raised  to  a  temperature  near  the  glass  transition 
temperature  Tg  of  the  membrane.  To  analyze  the  subsequent  development  of  stress  and  strain  in  the  two 
materials,  we  use  the  following  pair  of  simple  constitutive  relations  (which  follow,  for  an  axisymmetric  system, 
as  the  leading  order  constitutive  relations  of  each  theory  developed  in  [1]): 


Ei 


SflRi  —  S,  +  — — — j  (  €RR  +  Vi €©q  )  , 

1  V i 

E 

Seei  =  Si  +  —— — 2  (  C0©  +  VitRR  )  , 

1  vi 

where  S,  is  the  in-plane  residual  stress  in  material  i,  given  by 


(A.l) 

(A.2) 


5,  =  Sf  - 


Ei 

1- Vi 


CH  ( T  -  To) . 


(A.3) 


Here,  we  have  assumed  that  the  residual  stress  can  manifest  either  as  an  intrinsic  stress  S[,  associated 
typically  with  microstructural  changes  occurring  in  a  material,  or  as  a  thermal  stress  that* arises  if  the 
temperature  T  is  different  from  the  temperature  T0  at  which  the  stress  and  strain  of  the  material  is  purely 
mechanical.  For  thermal  stress,  a,  is  the  coefficient  of  thermal  expansion  (CTE)  of  material  i. 

We  take  the  membrane/mandrel  laminate  configuration  at  the  elevated  temperature  Tg  as  the  strain-free 
reference  configuration,  and  compute  subsequent  deformations  from  this  state.  We  assume  there  to  be  no 
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intrinsic  stresses  in  either  of  these  materials,  so  that  after  cooling  to  ambient  temperature  T  the  residual 
stresses  in  each  are  purely  thermal,  given  by 


S°  =  — 


Ea 


-as  ( T-Tg ),  and 


q°  —  Em 

bm  1-14 


am  (T  Tg). 


(A-4) 


Here,  the  s-subscript  refers  to  the  membrane  (since  it  will  be  the  substrate  for  an  optical  coating),  and  the 
m-subscript  refers  to  the  mandrel.  The  o-superscript  distinguishes  variables  prior  to  coating  and  releasing 
the  coated  membrane  from  the  mandrel. 

The  mandrel  is  assumed  to  undergo  a  free  contraction  during  cool-down  to  a  stress-free  state  at  the 
ambient  temperature  T,  so  that  for  i  =  m  the  left-hand  sides  of  equations  (A.l)  and  (A.2)  are  zero,  i.e., 

^RRm  =  Seem  =  hence 

0  =  ,  ™2  ( eRR  +  Um€e& ) »  (A- 5) 

0  =  Sm  +  :  (  £00  +  Vmt°RR  )  •  (A.6) 

These  equations  are  easily  solved  for  the  two  strain  components,  yielding 

e°RR  =  £©e  =  -!&  =  am  {T  ~  Ts)  >  .  (A'7) 

where  the  final  result  made  use  of  the  second  equation  of  (A.4).  Substituting  these  results  for  the  strain 
components  into  (A.l)  and  (A.2),  with  i  =  s,  then  yields  the  stress  components  in  the  membrane  substrate, 
i.e., 


S°rrs  =  S°ees  =  S°  +  (1  +  v.)  am  (T  -  Tg) 


Es  as  (T  —  Tg)  +  7—77  (T  —  Tg) 


1-^5 

Es 


1  -V8 

(as  —  am)  {T  —  Tg) , 


(A.8) 


1 

where  the  first  equation  of  (A.4)  was  used  to  get  the  second  line.  For  convenience,  we  write  the  last  equation 
as 

(A-9) 


Es 

Ss  = 


where 

and 


1  -  vs 


—  QO  _  CO 

S  —  dRRs  —  De&si 


es  —  ( cts  otm )  (T  Tg) 


(A.10) 
(A. 11) 


is  the  recoverable  strain,  i.e.,  the  strain  the  membrane  would  experience  if  released  from  the  mold  and  allowed 
to  come  to  a  stress-free  state. 

As  an  example,  we  consider  a  CP1-DE  membrane  material  (manufactured  by  SRS  Technologies,  Inc, 
Huntsville,  AL)  cast  on  an  aluminum  mandrel.  SRS  has  provided  the  following  material  properties  for 
CP1-DE: 


Es  =  2.172  GPa,  (A.12) 

v,  =  0.34,  (A. 13) 

as  =  51.2  ppm/  °C,  (A.14) 

ps  =  1.434  g/cc,  (A.  15) 

Tg  =  263  °C.  (A. 16) 
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From  the  Mat  Web  website,  we  find  that  aluminum  alloys  have  CTE’s  of  around  24  ppm/  °C.  Assuming  the 

final  temperature  to  be  an  ambient  room  temperature  of  T  =  20  °C,  the  recoverable  membrane  strain  from 
(A. 11)  is 

e*  =  (qs  -  Qm)  AT  =  27.2  x  10~6/°C  x  (-243  °C)  =  -0.00661,  (A.17) 

or  roughly  0.7%  strain.  From  (A.17)  and  (A.9)  the  residual  thermal  stress  in  the  membrane  is  then 
E 

S‘  =  ~\~^a  ts  =  3-291  X  109  x  000661  =  21.754  x  106  Pa  =  21.754  MPa,  (A.18) 

which  we  note  is  positive,  hence  a  tensile  stress. 

Suppose,  now,  that  a  single  optical  coating  is  applied  to  the  membrane,  and  assumed  to  be  perfectly 
bonded  to  the  membrane  (which  itself  remains  bonded  to  the  mandrel).  For  generality,  we  assume  the 
coating  is  placed  on  the  membrane  at  a  temperature  Tc,  either  above  or  below  the  temperature  T  at  which 
the  mandrel  is  stress-free.  When  the  coated  membrane/mandrel  is  returned  to  the  temperature  T,  the  strain 
components  are  given  by  (A.7),  assuming  the  mandrel  to  again  be  stress-free,  but  with  Tg  replaced  by  T 
Substituting  these  strain  components  in  (A.l)  and  (A.2),  now  with  i  =  c,  yields  the  stress  components  in 
the  coating  at  temperature  T: 


SRRc  -  Seec  =  Si  ~  (<*c  -  am)  (T  -  Tc) ,  (A.19) 

where  5/  is  the  intrinsic  stress  in  the  coating. 

When  the  coated  membrane  is  released  from  the  mandrel,  equations  (A.l)  and  (A.2)  again  apply  but 
with  the  residual  stresses  given  by  (A.8)  for  the  membrane,  and  by  (A.19)  for  the  coating,  i.e., 

Ss  =  s°RRs  =  5©0S  =  “• Qm)  (r  “  T9 ) .  (A. 20) 

Sc  =  S%Rc  =  5©0C  =  Src  -  - (ac  -  am)  (T  -  Tc) .  (A.21) 

If  the  coated  membrane  is  first  attached  to  a  rigid  boundary  ring,  and  then  released  from  the  mandrel, 
the  boundary  conditions  it  must  satisfy  are  those  of  a  rigidly  clamped  edge.  From  the  geometrically  linear 
theory  of  a  stress-coated  net-shape  film  with  a  clamped  boundary,  the  “on-design”  residual  coating  stress 
maintaining  the  initial  parabolic  shape  of  the  coated  membrane  after  releasing  it  from  the  mandrel  is  derivable 
from  equation  (4.171),  i.e., 

Af  -  2 aF*  (p  +  70  g)  =  0,  (A.22) 

where  we  have  used  (4.64)  and  (1.16)  to  replace  k0  in  terms  of  the  /-number.  For  a  two-layer  system  with 
no  pressure  difference,  and  ignoring  the  effects  of  gravity,  this  reduces  to 


Af  —  hcSc  +  hsSs  —  0,  (A. 23) 

where  hs  is  the  membrane  thickness  and  hc  is  the  coating  thickness,  hence  the  required  coating  stress  is 


Sc  = 


(A.24) 


Assuming  for  the  present  that  the  membrane  was  coated  at  the  temperature  Tc  =  T  where  the  mandrel  was 
stress-free,  the  thermal  term  in  (A.21)  is  zero,  hence  the  residual  coating  stress  is  purely  intrinsic: 


(A. 25) 


Thus,  for  a  1  pm  thick  coating  on  a  10  pm  thick  membrane,  and  a  CTE  mismatch  stress  in  the  membrane 
given  by  (A.18),  the  required  coating  stress  would  be 


10 

x  21.754  MPa  =  -218  MPa. 


(A.26) 
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Suppose,  however,  that  this  prescription  stress  level  Si  is  not  precisely  met,  but  acquires  some  value 
Sc  —  (1  +  cj>)S [  =  -(14-  4>)hsSs/hc  where  (j>  is  some  number  satisfying,  say,  -0.1  <<j><  0.1.  Then  we  have 

AT  =  -<j>hsSs  ±  0.  (A.27) 

Under  these  circumstances  it  may  be  feasible  to  apply  a  pressure  difference  to  null  out  the  displacement,  in 
which  case  the  condition  (A.22)  for  no  deformation  becomes 

A r  -  2 aF*p  =  0,  (A.28) 


neglecting  gravity.  According  to  (A.27)  and  (A.28),  applying  a  pressure  difference  of 

_  Af  _  ,  hsSs 

P  ~  2aF*  ~  ~(p2aF # 


(A.29) 


will  thus  maintain  the  initial  parabolic  shape.  For  a  one-meter  diameter  f/2  mirror  with  membrane  CTE 
mismatch  stress  21.754  MPa,  as  in  (A.18),  and  undercompensated  by  a  coating  stress  that  is  90%  of  the 
prescribed  value  (<f>  =  -0.1),  the  required  pressure  difference  would  be  very  small:  p  =  10.9  Pa,  or  about 
0.0016  psi. 
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